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Abstract The physics of collisionless shocks is a very broad topic which has been studied 
for more than five decades. However, there are a number of important issues which re- 
main unresolved. The energy repartition amongst particle populations in quasiperpendicular 
shocks is a multi-scale process related to the spatial and temporal structure of the electro- 
magnetic fields within the shock layer. The most important processes take place in the close 
vicinity of the major magnetic transition or ramp region. The distribution of electromag- 
netic fields in this region determines the characteristics of ion reflection and thus defines 
the conditions for ion heating and energy dissipation for supercritical shocks and also the 
region where an important part of electron heating takes place. In other words, the ramp 
region determines the main characteristics of energy repartition. All of these processes are 
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crucially dependent upon the characteristic spatial scales of the ramp and foot region pro- 
vided that the shock is stationary. The process of shock formation consists of the steepening 
of a large amplitude nonlinear wave. At some point in its evolution the steepening is arrested 
by processes occurring within the shock transition. From the earliest studies of collisionless 
shocks these processes were identified as nonlinearity, dissipation, and dispersion. Their 
relative role determines the scales of electric and magnetic fields, and so control the char- 
acteristics of processes such as of ion reflection, electron heating and particle acceleration. 
The determination of the scales of the electric and magnetic field is one of the key issues in 
the physics of collisionless shocks. Moreover, it is well known that under certain conditions 
shocks manifest a nonstationary dynamic behaviour called reformation. It was suggested 
that the transition from stationary to nonstationary quasiperiodic dynamics is related to gra- 
dients, e.g. scales of the ramp region and its associated whistler waves that form a precursor 
wave train. This implies that the ramp region should be considered as the source of these 
waves. All these questions have been studied making use observations from the Cluster 
satellites. The Cluster project continues to provide a unique viewpoint from which to study 
the scales of shocks. During is lifetime the inter-satellite distance between the Cluster satel- 
lites has varied from 100 km to 10000 km allowing scientists to use the data best adapted 
for the given scientific objective. 

The purpose of this review is to address a subset of unresolved problems in collisionless 
shock physics from experimental point of view making use multi-point observations onboard 
Cluster satellites. The problems we address are determination of scales of fields and of a 
scale of electron heating, identification of energy source of precursor wave train, an estimate 
of the role of anomalous resistivity in energy dissipation process by means of measuring 
short scale wave fields, and direct observation of reformation process during one single 
shock front crossing. 

Keywords Collisionless shocks • waves in plasmas ■ nonstationarity ■ shock scales ■ plasma 
heating and acceleration ■ wave-particle interactions 
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1 Introduction 

Collisionless shocks (CS) are ubiquitous in the universe. They play an important role in the 
interaction of the solar wind with the planets ( Russell 1 985 , 1 977 ; Greenstad t and Fredricks| 
|1979; Ness et al.|1974| [T~981 ), they also are supposed to have vital role in fundamental as- 
trophysical problems such as cosmic ray acceleration (Krymskii 1977[|Axford et al.|1977| 




Dynamic quasiperpendicular shock 



c 
o 



1 - 



non-relativistic 



£ I0" 3 i 

CD 
C 
CD 



Solar 



10 



-9- 



SNRs 
/irial 



? 



10 



-4 



10 



-2 



102 y S h(3si 



Fig. 1 (courtesy of A. Spitkovsky). Parametric range of observations of collisionless shocks associated with 
different astrophysical objects. 



Bell 1978a b, Blandford and Ostriker 19781. CS's are of crucial importance for understand- 
ing physical processes in the vicinity of such astrophysical objects as supernova remnants 
(Koyama et al. 1995; Ba mba et al.|2 003 1, plasma jets ( |Piran|2005] >, binary systems and or- 
dinary stars. In spite of this great variety of CS in the Universe only those shocks in the 
Solar system can be probed using in-situ observations. Moreover, comprehensive in-situ 
data exist only for interplanetary shocks and planetary bow shocks, however, it is worth not- 
ing that some astrophysical shocks are similar to those in the solar system. As was noted 
by [Kennel et aL|( |1985| 'The density, temperature and magnetic field in the hot interstellar 
medium are similar to those in the solar wind, and the Mach numbers of supernova shocks 
at the phase when they accelerate the most cosmic rays are similar to those of solar wind 
shocks'. Astrophysical shocks associated with different objects exhibit large differences in 
the parameters that characterise them. Figure [T] shows the variation of astrophysical shocks 
as a function of magnetisation (Y-axis), determined as 1/Ma where Ma is the Alfven Mach 
number (the ratio of the upstream flow velocity to the characteristic velocity of propaga- 
tion of magnetic perturbations in a plasma or Alfven velocity) and the characteristic plasma 
pressure to magnetic pressure ratio (X-axis) where 7^ is the ratio of the upstream flow ve- 
locity to the velocity of light and /3 S /, the ratio of total plasma particle thermal pressure to 
the magnetic field pressure in the reference frame of the upstream flow. Collisionless shocks 
associated with different astrophysical objects such as Supernovae Remnants (SNR), Active 
Galactic Nuclei jets (AGN), Pulsar Wind Nebulae (PWN), and Gamma Ray Bursts (GBR) 
are indicated. 

As can be seen from this figure, different 'families' of shocks occupy different regions 
of parameter space. One can also see that the parameters of SNR shocks are quite similar 
to those of solar system shocks. It allows one to suggest that the studies of solar system 
shocks, and in particular Earth's bow shock, represents an interest for wider scientific com- 
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munity than only for geophysics. The majority of astrophysical and Solar system shocks are 
developed in magnetised plasmas. 

Collisional shocks have been studied for many decades, beginning with the earliest ob- 
servations of Mach (Mach and Arbes 1886; Mach and Salcher 1887 ). A shock occurs when 
an obstacle finds itself immersed in a supersonic gas flow. Before reaching the surface of the 
obstacle the flow should be decelerated to velocities lower than the velocity of sound so that 
it may flow around the body. This process of flow deceleration and the redistribution of its 
directed energy occurs over distances of the order of the collisional mean free path of gas 
particles and the energy difference. The energy taken from the flow during deceleration is 
mainly transformed into the thermal energy of the gas as it is heated. As a result the sound 
velocity in the gas increases and, after the shock transition, becomes larger than the remain- 
ing directed velocity of the flow so that the motion downstream of the shock is subsonic. 
Thus the shock represents the transition from supersonic directed motion to subsonic in the 
reference frame of the obstacle immersed into the flow. 

The notion of the collisionless shock was introduced by several authors in the late 50's 
( |Adlam and Allen|[l958l |Davis et al.|1958|[Sagdeyev|[T965i >. The modern form of the de- 
scription was presented in an almost complete form in the famous review paper by Sa gdeev| 
( 1966). The first problem to be overcome is related to the the existence of a shock. For col- 
lisional shocks (as mentioned above) the shock thickness is related to the mean free path of 
the gas particles. However, in space plasmas the mean free path can be as large as 5AU (i.e. 
similar to the distance of Jupiter from the Sun)! So, how can a shock exist whose width is 
much smaller than the mean free path? Historically, a very similar problem first appeared in 
laboratory devices and only later in space plasmas. However the crucial issue in both cases 
is ex actly the same (|Paul et al.||1965| |Kurtmullaev et al.||1965| |1967| |Ascoli-Bartoli et al.| 
|1966[|Goldenbaum and Hintz|1965fr . 

The solution proposed initially relied on the processes of anomalous dissipation, namely, 
anomalous resistivity. Thin shock transition contains quite strong variations of the magnetic 
field components perpendicular to its normal. This implies that there is a very intense cur- 
rent inside the shock transition layer. The current carriers, charged particles move with re- 
spect to another. The plasma state supporting these intense currents is, in general, unstable. 
The instabilities in the plasma result in the generation of intense waves. Wave generation 
opens new channel of impulse and energy exchange between the different populations of 
plasma particles. For instance, current carrying electrons can emit/generate the waves, and 
these waves can be absorbed by ions. This generation-absorption exchange using waves as 
a transmission media between plasma components leads to an exchange of energy and pulse 
between them. Typically, the characteristic scale of energy exchange between different par- 
ticle populations can be much smaller than the mean free path of particles. As a result the 
characteristic scale of the dissipation process can be determined by this anomalous dissipa- 
tion. Thus the principal difference between collisional and colisionless shocks is the change 
of the dissipation scale that is determined by additional process involved, but the nature of 
the transition and its characteristics remain very similar. In both cases the shock redistributes 
the directed bulk plasma flow energy to plasma thermal energy. However, the dissipation rate 
and characteristic scales of collisionless shocks are determined by the anomalous process of 
energy dissipation. The notion of anomalous resistivity was already well known and widely 
used in plasma physics. The theory of anomalous resistivity based on current instabilities 
and the generation of ion-sound waves was directly applied to the theory of collisionless 
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( 1976| l. Later this idea was further developed in series of papers by [Pa 



3 a I, who noticed that in the case of currents perpendicular to the magnetic 
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propagate almost perpendicular to the magnetic field and the theory of anomalous resistivity 
in this case should account for these rather than ion-sound waves. 

From the earliest experimental studies of shocks in space and laboratory plasmas it was 
found that the characteristics of the shocks observed can be quite different even in the range 
of parameters that correspond to solar system shocks and those in laboratory plasmas. There 
were observations of quite small scales for the ramp with much longer precursor wave train, 
there were shocks consisting of a long transition region with large amplitude structures in 
the magnetic field filling a very large area in space. These early observations gave rise to 
attempts to classify shocks. 

The first systematic classification was proposed by Formisano (1985 ). He noticed that 
there are three basic parameters of the upstream flow that are important for the classifi- 
cation. These are the angle between the magnetic field and shock front normal, 6g n , the 
plasma beta /3 i.e. the ratio of total particle thermal pressure to the magnetic field pressure 
P = 8-rmT/B 2 , where T is the total plasma temperature, n is the plasma density, and B the 
magnitude of the magnetic field; and the magnetosonic Mach number Mms = (Vup/Vms), 
where V up is the normal component of the velocity to the shock, Vm s is me velocity of the 
magnetosonic wave propagating in the same direction as the shock. Later this classification 
was slightly modified and is used in the form proposed by Kenne l et al.| ( |1985] l. This paper 
divides shocks in to two broad classes that are related to the ion dynamics, namely, quasi- 
parallel and quasiperpendicular. The characteristic feature of the first group is that the ions 
that are reflected from the shock front can freely propagate to the upstream region. These 
shocks correspond to the range of angles between the magnetic field and the normal vector 
to the shock front 8g n < 45°. In the second group, quasiperpendicular, part of the ion pop- 
ulation is reflected. After reflection they turn back to the upstream region and can gain extra 
energy due to the acceleration by inductive electric field tangential to the shock surface and 
perpendicular to the magnetic field. Then they can cross the shock front ( |Woods|1969l[l97T] 
Sckopkeet al.|1983| >. This process can occur when 8g n > 45°. 



Low Mach number shocks could dissipate the necessary energy entirely through some 
anomalous resistivity within the current-carrying shock layer. The right-hand fast magne- 
tosonic/whistler waves have phase and group velocities that increase with decreasing wave- 
length beyond the fluid regime. Thus, steepened fast mode shocks are expected to radiate 
short wavelength waves, and hence energy, into the unshocked oncoming flow. The shortest 
wavelength capable of standing in the flow then forms a "precursor wavetrain" that has been 
observed at these sub-critical shocks (Mellott 1985 ) and as we shall show later this pro- 
cess occurs in supercritical shocks also. However, above a critical Mach number, anomalous 
resistivity within the layer carrying the limited shock current is unable to convert the re- 
quired amount of energy from directed bulk flow into thermal energy. At quasi-perpendicular 
shocks, where the magnetic field in the unshocked region makes an angle 8g n > 45° with 
the shock propagation direction (the shock "normal" n), a fraction of the incident ions are 
reflected by the steep shock ramp as described above. They gyrate around the magnetic field 
and gain energy due to acceleration by the transverse motional electric field (—V x B). Re- 
turning to the shock layer they have sufficient energy to pass through into the downstream 
shocked region (Woods 1969, 1971; Sckopke et al.|1983] >. The separation of ions onto two 
groups, crossing front directly and after the reflection results in the dispersal of particles in 
velocity space. Group of reflected particles is separated from the bulk ion population due 
to an increase in peculiar velocity relative to the bulk motion. This process corresponds to 
the kinetic "heating" required by the shock jump conditions and it ensures the major part 
of energy dissipation necessary for directed energy transfer to thermal energy of plasma ion 
population. The process of ion thermalization takes place on rather large scale downstream 
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of the shock front. The spatial length of the transition to ion thermal equilibrium can be 
treated in a similar fashion to that of the shock front thickness in collisional shocks. De- 
tailed measurements of ion distributions onboard ISEE mission resulted in establishing all 
major characteristics of this process ( Sckopke et al. 1983 ). This result is probably one of the 
most important obtained in this outstanding program. In theory this critical Mach number 
corresponds to the multi-fluid hydrodynamic limit in which the resistivity and viscosity can 
nolonger provide sufficient dissipation l |Coroniti|1970| >. The reflection occurs on sufficiently 
smaller scales than thermalization due to a combination of magnetic forces and an electro- 
static cross-shock potential. The main potential, which corresponds to the frame-invariant 
E • B electric field, is known as the deHoffmann-Teller potential ( [de Hoffmann and Teller| 
1950; Goodrich and Scudder 1984 ). It results directly from the leading electron pressure 
gradient term in the Generalised Ohm's Law (Scudder et al. 1986b). However in more de- 
tailed two-fluid descriptions, the quasiperpendicular shock has fine structure that depends 
upon the characteristics of the nonlinear shock profile (Galeev et al. 1988, 1989; Gedalin 
\TWf\ [Kras noselskik h et al.||2002] ). In this paper we shall present the results of the stud- 
ies of quasiperpendicular shocks only. Quasiparallel shocks will be discussed in a paper by 
Burgess and Scholer (this review) . 

The basic ideas of the shock formation can be understood by considering the propa- 
gation and evolution of large amplitude wave. In gas dynamics the wave corresponds to a 
sound wave whose evolution in terms of gas dynamics leads to the formation of disconti- 
nuities. In reality, however, narrow transition regions are formed in which the dynamics are 
dominated by dissipative processes. In plasmas, the characteristics of the main shock transi- 
tion are determined not only by an interplay between nonlinearity and dissipation, but also 
by another important physical effect, the wave dispersion. It is well known that a subcriti- 
cal shock has a nonlinear whistler wave train upstream of its front (Sagdeev 1966; Mellott 
|1985| l. The presence of whistler/fast magnetosonic precursor wave trains in supercritical 
shocks as well, was experimentally established in Krasnoselski kh et al.| (1991) ; |Bali khin 
|et al.| ( |1997b| l; |Oka et al] ( |2006[ >. These whistler waves have rather large amplitudes and 
their role in energy transformation and redistribution between different particle populations 
and in the formation of the structure of the shock front is still an open question. The major 
transition of such a dispersive shock, the ramp, may behave in a similar fashion to either 
the largest peak of the whistler precursor wav e packet ijK arpman e t al.|1973j|Sagdeevf l966 ; 
|Kennel et al.|1985[[Gaieev et al.|1988"l[T989l [Krasnoselskikh et al.|2002| l or the dissipative 
shock region in which the major dissipation due to current driven instabilities occurs. The 
nonlinear steepening process can be described as the transfer of energy to smaller scales. 
The steepening can be terminated either by collisionless dissipation, as described above, or 
by wave dispersion. Typically the dissipative scale L d exceeds the dispersive one L(H S p, the 
former is reached first and further steepening can be prevented by the dissipation that takes 
away energy. When steepening is balanced by dissipation, a dissipative subcritical shock 
forms. Most subcritical collisionless shocks observed in situ are supposed to be dissipative 
even though dispersive processes play a role in forming a dispersive precursor wave train. 
Such a shock is characterized by a monotonic transition in the magnetic field (magnetic 
ramp) of the width L d . The dissipative length is determined by the most important anoma- 
lous dissipative process. Its major features are the generation of intense short-scale waves 
and their dissipation. This form of the evolution of a nonlinear wave takes place at low 
Mach numbers. However, if the nonlinearity is strong enough (as determined by the veloc- 
ity and density of the incoming flow), dissipation is not capable of stopping the steepening, 
and the gradients continue to grow then energy transfer to smaller scales continues and the 
characteristic scale of the transition can become smaller. The next process that comes into 
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play to counterbalance the steepening is dispersion. Dispersion becomes important when the 
gradients become comparable with the dispersive scale L disp . In this case the shock front 
structure becomes multi-scale. The steepening is prevented by short-scale dispersive waves 
which are able to propagate away from the evolving shock front. These waves effectively 
remove some part of the energy and, most importantly, restrain further growth of the gra- 
dients. For perpendicular shocks the phase velocity of the dispersive waves decreases with 
decreasing scale and a wave train is formed downstream of the magnetic ramp. For shocks 
with a more oblique geometry (quasiperpendicular shocks) the phase velocity increases with 
decreasing spatial scale and an upstream wave train is formed. The upstream wave precur- 
sor is approximately phase standing in the upstream flow. Its amplitude decreases with the 
distance from the shock ramp due to dissipation processes as was discussed in the early 
theoretical papers describing subcritical shocks (Sagdeev 1965a, 1966). 

The transition to reflection shock takes place when downstream bulk velocity reaches 
the downstream ion-sound speed. Supercritical reflection shocks have a more complex struc- 
ture in comparison to subcritical shocks. In quasi-perpendicular shocks the upstream mag- 
netic field does not allow reflected ions to travel far upstream before turning them back to 
the shock front. The upstream region in which the beam of reflected ions perform part of 
their Larmor orbit before being turned back to the shock is called foot. The foot region is 
characterised by a 15-20% increase in the magnitude of the magnetic field. The consider- 
ation of a Larmor orbit of a reflected ion gives relatively accurate estimate of the spatial 
size of the foot Lf = 0.68-R/^ smdg n where Rn is the gyroradius of ions moving with 
the velocity equal to normal component of the velocity of upstream flow (Woods 1969; 
|Livesey et al.|19 84). The coefficient 0.68 corresponds the case of to specular reflection. For 
non-specular reflection this relation is slightly modified ( |Gedalin|~1 996). Downstream of the 
quasi-perpendicular shock's main transition the joint gyration of the bulk plasma ions and 
beam of reflected ions leads to an overshoot-undershoot structure. Again, the size of this 
overshoot/undershoot can be estimated in a straightforward manner in terms of the ion gy- 
roradius. However, the main transition layer lies between the foot and the overshoot. This 
is the region where the most dramatic changes in the plasma parameters occur. In a super- 
critical, quasi-perpendicular shock this layer is characterised by the steepest increase of the 
magnetic field referred to as the ramp. The change of the electrostatic potential, reflection of 
ions, and electron thermalisation take place within the ramp and its spatial scale determines 
the major physical processes within the shock and the mechanisms for the interaction of the 
shock front with the incoming electrons and ions. For instance, several theoretical models 
suggest that in the ramp of high Mach number shocks very small scale electric fields can 
be present (Krasnoselskikh 1985; Galeev et al. 1988, 1989). There are several critical issues 
regarding supercritical quasiperpendicular shock physics for which alternative explanations 
for the observational features of the shock front have been proposed. Theoretical considera- 
tions ( jGaleev et al.|1 988, 1989; Krasnoselskikh et al. 2002), that treat the supercritical shock 
front as being similar to a nonlinear dispersive wave, predict that the ramp scales (gradients) 
should decrease with increasing Mach number, eventually reaching characteristic values as 
small as several electron inertial lengths L e = c/lj p . Moreover, after some critical Mach 
number corresponding to nonlinear whistler critical Mach number whose value is approx- 
imately 1 .4 times the linear whistler critical Mach number, the shock should become non- 
stationary. These critical Mach numbers determine the characteristic flow velocities when 
they become larger than the maximum velocity of a linear or nonlinear whistler wave prop- 
agating upstream along the shock front normal. Many computer simulations ( Schol e~and| 
Matsukiyo 2004; Matsukiyo and Scholer 2006) come to the conclusion that the thickness of 
the ramp is determined by the dissipative process due to either the modified Buneman in- 
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stability (MBI) or the modified two stream instability (MTSI). Both theoretical studies and 
computer simulations have associated pitfalls. Theoretical models can not accurately take 
into consideration the presence of reflected ions whilst simulations are carried out with an 
unrealistic ratio of the plasma frequency to gyrofrequency that strongly changes the ratio 
of electric to magnetic wave fields and often with an unrealistic ion to electron mass ratio. 
Both theoretical models and simulations predict the transition to nonstationary dynamics. 
However, they strongly differ in the determination of the scales of the electric and mag- 
netic fields in the ramp region, in the energy sources for the upstream whistler waves that 
form the precursor wave train, and in the characteristics of the shock dynamics when it be- 
comes nonstationary. For these reasons experimental studies of these questions are crucial 
for our understanding of the physical processes in quasiperpendicular collisionless shocks. 
Our Review aims to report the studies of all these questions making use mainly of Cluster 
measurements (adding some other data where it is necessary, in particular THEMIS data in 
studies of magnetic field scales of shocks). 

The first critical issue we shall address is magnetic ramp width and spatial scale. The 
main motivation for the study of the magnetic ramp width L r is that it is this scale that deter- 
mines the nature of the shock, i.e., the dominant physical processes that counteract nonlinear 
steepening. The shock width can be determined either by the solitary structure of nonlinear 
whistler slightly modified by the presence of reflected ions or by characteristic anomalous 
resistivity scale associated with one of instabilities mentioned above. For instance if it is in- 
deed the case that the ramp width increases with increasing Mach number as concluded by 
Bale et al. (20031, then the evolution of nonlinear whistler waves must be excluded from the 
processes that are involved in the formation of supercritical shocks. The characteristics of 
the major transition within the shock in which the flow deceleration and the magnetic field 
and electrostatic potential variations take place are determined by the interplay between non- 
linearity, dispersion and dissipation. The presence of a population of reflected ions makes it 
difficult to construct a reliable theoretical model based on an analytical or semi-analytical 
description. However, the establishment of the scales of this transition allows one to deter- 
mine the characteristics of the dominant physical processes in play. The ramp thickness is 
also crucial for a redistribution of energy between electrons and ions. An important char- 
acteristic involved in this process is the gradient scale of the transition. Two reasons cause 
the need in introducing this ramp gradient spatial scale. The first is the interaction between 
the incoming electrons and the electromagnetic field at the shock front. As was shown by 
Balikhin et aT] (19931 ; [Balikhin and Gedalin| ( fl994l i; |Gedalin et~aT1 l | 1995 alb) , and |Balikhin1 



et al. ( 1998[ l, an important effect of this interaction is the possible violation of adiabaticity 



even in the case when the width of the magnetic ramp considerably exceeds the formally 
calculated electron gyroradius. Two very different scenario of electron heating can occur 
depending upon if the conditions for adiabaticity are satisfied or violated. This effect is cru- 
cially dependent upon the ramp spatial scale. The change between adiabatic/non-adiabatic 
regimes is related to the ability of the nonuniform electric field within the ramp to rectify 
the electron motion and increase their effective gyration radius ijBalikhin et al.|1993"j |Ba- 



|likhin and Gedalin|1994HGedalin et al.|1995a|b[|Balikhin et al.|1998| >. The parameter that 



determines the transition from adiabatic to nonadiabatic motion of the electrons is the in- 
homogeneity of the magnetic field. The characteristic spatial scale of such a layer may be 
defined as the product of the change in magnetic field AB normalised to the upstream field 
Bo and the spatial distance dx over which this change occurs i.e. l gr = Lb t {Bo/ AB). 
To illustrate the effect of this parameter, one can consider two cases. For a weak shock for 
which (B^/Bq) ^ 1.2, here B^ is the magnetic field magnitude after the shock transition 
(downstream), and whose ramp width is of the order 5-6 Rlc ( electron Larmor radii) the 
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electron motion will be adiabatic. However, in a stronger shock of similar magnetic ramp 
width for which the maximum magnetic field observed in the overshoot exceeds that of the 
upstream field by a factor 5-6 the electron behaviour becomes non-adiabatic. This makes it 
necessary to carry out the statistical study of the magnetic ramp spatial scale in addition to 
the ramp width (size). The ramp width and its gradient scale are also important for the prob- 
lem of stability of the ramp region of the shock front. According to Krasnoselskikh ( 1985 1, 
|Galeev et al] ( |1988| |1989[ l and |Krasnoselskikh et al.| l |2002| l the nonlinear whistler wave 
structure becomes unstable when the characteristic gradient exceeds some critical value. It 
was suggested by Krasnoselskikh ( 1985 ) that it takes place when the Mach number becomes 
equal to nonlinear critical whistler Mach number. When this happens dispersive process can 
no longer counterbalance the nonlinearity and the shock front overturns. Thus, the charac- 
teristic gradient scale provides a rather universal characteristic of the degree of steepness 
of the shock front. Thus its determination completed a comprehensive statistical study of 
the magnetic ramp spatial gradient scale in addition to the ramp width (size). Many papers 
were devoted to the magnetic field structure of collisionless shocks (e. g.|Russell et al.|1983| 
Krasnoselskikh et al.|[T99Tj |Farris e"TaT1[T99T) |Newbury and Russell||1996| |Hobara et al'| 



2010| Mazelle et al.|2010 1. In particular, the spatial scales of its various regions have been 
comprehensively investigated (Bal ikhin et al.|1995||Farris et al.|1993[|Hobara et al.|2010| 
Maz elle et al.|2 010). ISEE and AMPTE measurements of the magnetic field profiles of the 
shock front structure led to evaluation of the scale sizes of the foot and overshoot regions 
that were supposed to be of the order of ion inertial length c/uj pi and 3c/ui pi respectively, 
here cu pi is the ion plasma frequency. The ramp scale has been estimated to be less than an 
ion inertial length with reports of one or two very particular shocks whose ramp scale was 



sufficiently smaller, of the order 0.lc/uj p i I Newbury and Russell 1996; Walk er et al.| 1999b} . 



We report here the statistical studies of scales based on papers by |Hobara et al.| l |2010| > and 
Maz elle et al.| ( (20T0| l. Another critical issue we shall address hereafter is the electric field 
distribution inside the ramp region. The energy transfer to smaller scales due to steepening 
can achieve the scales comparable to electron inertial scale where the whistler waves be- 
come quasi-electrostatic. In nonlinear-dispersive scenario of the shock front description the 
field can have multiple short scale electric field spikes. The experimental study can answer 
the question do they exist or not and to determine their characteristics. Studies of the mag- 
netic field profile across the terrestrial bow shock significantly outnumber those based on 
electric field measurements. Despite the fundamental effect that the electric field has on the 
plasma dynamics across collisionless shocks, the complexity of the interpretation of elec- 
tric field data has impeded studies of the electric field structure within the shock front. It 
is worth noting that only a handful of studies are dedicated to the electric field structure 
within the shock front (Heppn er et al.|1978)|Formisano|1982)|Scudder et al.|1986a[|Wy gant 
let al.| 19871 |Balikhin et al.|20021|Walker et al.|2004||Balikhin et al.|2005HBale and M ozer 
[20071 [Hobarae t al. 2008] |Dimmock et al.||2012||Bale et al.|2008| > In contrast, there have 
been very few reports regarding the scale lengths of features observed in the electric field 
at quasi-perpendicular shocks. The scale size over which the potential varies at the front 
of a quasi-perpendicular bow shock is an issue that requires resolving in order to gain a 
full understanding of the physical processes that are occurring in the front. Several different 
points of view have been published on the relationship between the scale size of the mag- 



netic ramp and that over which the change in potential is observed. Some studies (Eselevich 
et al.|1971t|Balikhin et al.TT993j [F ormisano and Torbert|1982[|Formisano|1982"l|1985t|Ba 



likhi n et al.|2002l |Krasnoselskikh|1985| |Leroy et al.| 19821 |Liewer et al.|1991| |Scholer et aL 



2003 ) have proposed that the spatial scale of electrostatic potential is of the same order or 
smaller than that of the magnetic ramp under certain conditions. Such shocks have been ob- 
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served in numerous experimental and numerical studies of quasi-perpendicular supercritical 
shocks. On the other hand Scudder (1995) claimed that the potential scale length is larger 
than that of the magnetic scale length. Actual measurements of the electric field variations 
within the bow shock are very sparse. The main reason for this is due to the difficulties en- 
countered when making electric field measurements. Only a small number of space-based 
measurements of the electric field during the shock front crossing have been reported. Hep- 
|pner et aL] ( |1978| ) reported observations of a short lived spikes in the electric field making 
use of ISEE measurements. However, being short duration, these features were not observed 
at every shock crossing. Subsequent investigations by Wygant et aT|l |1987| > have shown the 
existence of spike-like features in the electric field both at the shock ramp and in the region 
just upstream. From the study of spin averaged ISEE-1 data, Formisano (1982) determined 
that the increase in the observed electric field E intensity began just upstream of the mag- 
netic ramp and lasted longer than the ramp crossing itself. Whilst the E-field intensity in the 
regions upstream and downstream of the shock could be interpreted as due to the V x B 
motion of the plasma the enhancement observed during the shock crossing must be due to 
the processes occurring within the shock front itself. In laboratory experiments where the 
conditions are not exactly the same as in space plasmas Eselevich et al.| l [T97T] l reported that 
the major change in potential across the shock occurs within the magnetic ramp region. 

Using data of numerical simulations, Lemb ege et al.| (19991 analysed the scale size of 
both the magnetic ramp region Lg r and the scale of the major change in potential in- 
side and around the ramp. They concluded that the scale lengths were of the same order, 
i.e. Lbt ~ L^. This view is supported by the simulations of Scholer et al. (2003 i. The 
latter authors also show that during the shock reformation process, the main potential drop 
occurs over several ion inertial scales in the foot region and they noticed that the steepened 
magnetic ramp region also contributes a significant fraction of the change in total poten- 
tial over much smaller scales, typically 5-10 Debye lengths. Despite the simulation shocks 
parameters are still rather far from observations (see Section[3]for more details) the tenden- 
cies in majority of simulations are well pronounced and are similar to those in laboratory 
experiments. Hereafter we report the observations of electric field spikes observed onboard 
Cluster satellites first reported by Walker et al. (2004) and the statistical study of their char- 
acteristics. The third important problem of quasiperpendicular shock physics addressed in 
this Review is the problem of electron heating. By contrast to ion heating problem well 
advanced due to detailed studies onboard ISEE mission, the electron heating problem has 
remained controversial. The action of shock quasistatic electric and magnetic fields on the 
electron population (which can have thermal speeds far in excess of the shock speed) is to 
inflate and open up a hole in the phase space distribution by accelerating (decelerating) in- 
coming (escaping) electrons ( Scudder et al.|1986c||Feldman et al.|1982[ ). This inflation in 
itself is reversible thus it is not dissipation or heating if other processes would not be in- 
volved. Irreversibility may be imposed if additional scattering would take place infilling the 
hole. If adiabatic invariant of electrons is conserved while electrons cross the shock front it 
can not happen. One should conclude that some non-adiabatic process should occur inside 
the shock front. One of the possibilities can be related to Debye-scale electric fields l jBale| 
|et al.| 1998] ). Another possibility is to suggest that the phase space inflation is indeed accom- 
panied by instabilities which could scatter the electrons. Demagnetisation of the electrons 
due to the strong gradients in the electric field as it was mentioned above (Balikhin and 
Gedalin 1994) or nonlinear wave phenomena Krasnoselskik h et al.| ( |2002} combined with 
wave particle interaction can offer alternative scattering processes. Thus the partition of en- 
ergy between ions and electrons is a complex, self-consistent multi-scale interplay between 
electron heating, magnetic/electric field profile, shock potential, and ion reflection. This in- 
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terplay remains poorly understood despite 40 years of research. That research has included 
detailed case studies ( Scudder et al.|1986b)>, statistics of the inferred potential and electric 
field structures (Schwartz *et al.|1988[ [Walker et al.|2004[ l, theoretical studies l |Galeev etaL] 
T988l|Gedalin|1997HKrasnoselskikh et al.|2002|> and increasingly sophisticated numerical 
simulations (Lembege et al. 2004; Scholer and Burgess 20061. Direct measurements of the 
thickness of the shock transition layer combined with the rapid simultaneous measurements 
of the electron distribution function can allow solving this long standing opened problem 
in shock physics. If the electron heating can be attributed to kinetic instabilities, the shock 
thickness will be measured in ion inertial lengths (c/w p i) (Papadopoulos 1985b; Matsukiyo 



and Scholer 2006 ). If such instabilities prove ineffective, above a second critical Mach num 
ber the shock steepening is expected to be limited by whistler dispersion and/or be unstable 
to shock reformation (Krasnoselskikh et al. 2002). Recent studies of the shock thickness 
(Hob ara et al.||20i0) jMazelle e t al. 2010 } do show scales comparable to whistler wave- 
lengths. These contrasted an earlier study ( |Bale et al.|2003| > reporting scalings that matched 
the gyro-scales of reflected ions. To date, studies have relied on the high temporal cadence 
available from magnetic or electric field experiments. However, field profiles provide only 
indirect evidence of the shock dissipation scales. A recent study ( |Lefebvre et al.|2 007 ) used 
sub-populations of electrons to determine the electrostatic potential profile at one shock, 
suggesting that it rose in concert with the magnetic field. In the work reported here, first 
published in (Schw artz et al.|2011[ l, the electron distribution function major characteristics 
are measured at sufficient cadence to reveal directly for the first time the scale of the electron 
temperature profile. Many shock crossings by Cluster satellites take place on the flanks of 
the magnetosphere that creates quite favourable conditions for the studies of the relatively 
narrow shock transitions allowing one to have many measurements on small spatial scale. 
Hereafter we present unprecedently rapid measurements of electron distribution moments 
that allow to shed new light on electron heating problem and its scales. 

The fourth problem presented in the Review, is closely related to the problem of mag- 
netic and electric field scales, is determination of the source of waves forming upstream 
precursor wave train. The presence of whistler/fast magnetosonic precursor wave trains in 
supercritical shocks was experimentally established in Bali khin et al.| ( |T997a| l;|Krasnosel- 
|skikh et al.|(T99T] l; |Oka et al.| ( (2006] . These whistler waves have rather large amplitudes and 
their role in energy transformation and redistribution between different particle populations 
and in the formation of the shock front structure is still an open question. The energy source 
responsible for the generation of these waves is the subject of active debate in shock physics 
(seejGaleev et al.|[T98"8"l [1989] |Krasnoselskikh et al.||20"0"2"[ |Matsukiyo and Scholer||20"06{ 



|Comi§el et al.||2011| >. Often the precursor waves are almost phase-standing in the shock 
frame. However, if they are generated by the ramp region as the dispersive precursor their 
group velocity can still be greater than zero in the shock reference frame, which would allow 
energy flow in the form of Poynting flux to be emitted towards the upstream of the shock 
transition. On the other hand, if the waves are generated by instabilities related to reflected 
ions their energy flux will be directed from the upstream region towards the shock ramp. 
The goal of Section]?] is to address this problem, to present the direct measurement of the 
Poynting flux of the upstream whistler waves aiming to establish the direction of the Poynt- 
ing flux. There are two different points of view on this subject also. It has been suggested 
that the shock front structure of quasi-perpendicular supercritical shocks is formed in a way 
similar to that of subcritical shocks ( |Galeev et al.|1988|[T989l|Krasnoselskikh et al.|2002) . 
In such scenario the precursor wave train is a part of the shock front structure emitted by the 
ramp region upstream due to positive dispersion of whistler waves. The observed dynamic 
features of shocks have also been studied extensively using computer PIC- or hybrid sim- 
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ulations, often with focus on the precursor wave activity and reflected ions ( [Hellinger and| 
Mangeney 1997; Hellinger et al. 2007; Matsukiyo and Scholer 2006). From a kinetic view- 
point, however, it may be argued that the shock-reflected ions change the physical picture 
and that the principal scales, temporal and spatial, could be determined by the characteris- 
tics of the reflected ion population (Bale et al. 2003 ). Upstream waves can then be generated 
due to counterstreaming ions and electrons in the shock front region, forming unstable par- 
ticle distributions with respect to so me wave modes (jPa padopoulos 1985b, Hellinger et al.| 
|2007| |Scholer and Matsukiyo|2004| |Matsukiyo and Scholer|2006| l. While this is probably 
the case for some higher frequency waves, we present here an analysis that leads to the con- 
clusion that the source of the upstream low frequency whistler waves is indeed related to 
the presence of the nonlinear ramp transition, emitting smaller scale dispersive waves to- 
wards the upstream flow. The existence of phase-standing upstream whistler waves depends 
on the value of the upstream flow speed Mach number relative to the phase velocity. If the 
Mach number of the shock does not exceed the nonlinear whistler critical Mach number 
M w = V w ,max/VA = l/2y / mi/m e cos 0g n , where V w ,max represents the highest possi- 
ble velocity of nonlinear whistler wave, then phase-standing (nonlinear) whistler wave trains 
can exist upstream of the shock l |Galeev et al.|1988||1989| Krasnoselskikh et al.|2002[ ). The 
results we report here were first published in ( [Sundkvist et al.|2012 1. Similar results were 
reported making use the Time Domain Sampling instrument (TDS) onboard Wind satellite 
(Wil son et al.|2012| for three of four crossings of interplanetary shocks. In one case the po- 
larization of waves was found to be different from whistler wave. Unfortunately, one satellite 
measurements do not allow to establish unambiguously the reason of this anomaly, it could 
associated with some particular perturbation in the solar wind. 

The fifth problem intimately related to previous is the problem of nonstationary dy- 
namics of high Mach number shocks. Shock waves are usually considered to be nonlinear 
waves that cause irreversible changes of state of the media and from macroscopic point 
of view they are stationary (for a review, see, e.g., |Tidman and Krall {T97T| ). However, in 



the very beginning of the collisionless shock physics [Paul et al. ( 1967 1 hypothesized that 
high-Mach-number shocks can be nonstationary, and the first unambiguous evidence of the 
nonstationarity was obtained by |Morse et al.| l fl972| > in laboratory experiments. New evi- 
dence of shock front nonstationarity was found in the 1980s. In particular, Vaisberg et al. 
( 1984) reported low frequency oscillations of the ion flux in the Earth's bow shock. Later 
jeg enal et al.|(1987) observed a similar phenomenon in the Uranian bow shock. In the very 
beginning of computer simulations of the collisionless shocks Biska mp and Welter] ( |1972| ) 
have observed the process of shock dynamic behaviour. The inflowing ions formed vortices 
in the phase space and dynamics of the front was definitely nonstationary. Later, numerical 
simulations performed by |Leroy et al.| lfl982) using 1-D hybrid code showed that the front 
structure of perpendicular shocks varies with time, for instance, the maximum value of the 
magnetic field exhibits temporal variations with a characteristic time of the order of ion gy- 
roperiod, the magnitude of these variations being about of 20% if the parameters are typical 
for the Earth bow shock (Ma = 8 and f3 e j = 0.6, where Ma is the Alfven Mach number, 
f3 C) i is the ratio of the thermal electron/ion and magnetic pressures). They also found that for 
Ma = 10 and /3 e i = 0.1 the ion reflection was bursty, oscillating between and 70-75%. 
Hybrid simulation of perpendicular shocks with very high Mach numbers carried out for the 
first time by Quest ( 1986) have shown that the ion reflection in the shocks can be periodic, 
the stages with 100% ion reflection alternating with the stages of 100% ion transmission. As 
a result, instead of a stationary structure, he observed a periodic wave breaking and shock 
front reformation. Later jHellinger et al. 1(2002") reexamined the properties of perpendicu- 
lar shocks with the use of the 1-D hybrid code and observed the front reformation for a 
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wide range of parameters if upstream protons are cold and/or Mach number is high. |Scholer| 
[etaT1 ( (2003] l and |Scholer and Matsukiyo] l |2004| ; [Matsukiyo and Scholer| ( |200o} in their 1-D 
full-particle simulations with the physical ion to electron mass ratio reproduced the refor- 
mation of exactly and approximately perpendicular high-Mach-number shocks in plasmas 
with /3i = 0.4 and demonstrated an importance of modified two-stream instability for the 
reformation process. Krasnoselskikh~| < p~985[ > and |Galeev et al.| ( p"98 8 19891 proposed mod- 
els describing the shock front instability due to domination of nonlinearity over dispersion 
and dissipation. This instability results in a gradient catastrophe within a finite time inter- 
val. Several aspects of the model, including the role of nonlinear whistler oscillations and 
existence of a critical Mach number above which a nonstationarity appears, were developed 
in further detail and more rigorously by Krasnoselski kh et al.| ([2002 ) and complemented 
by numerical simulations with the use of the 1-D full particle electromagnetic code with a 
relatively small ratio of electron and ion masses, m e /m,i = 0.005. It was also shown that 
the transition to nonstationarity is always accompanied by disappearance and re-appearance 
of the phase-standing whistler wave train upstream of the shock front. Moreover, for large 
Mach numbers the nonstationarity manifests itself as a periodic ramp reformation, which 
influences considerably the ion reflection, in particular, the reflection becomes bursty and 
sometimes the ions are reflected from both old and new ramps simultaneously. The four- 
spacecraft Cluster mission gave much more new opportunities for experimental studies of 
the shocks. The first examples of some aspects of shock nonstationarity were presented by 
|Horbury et al.| ( [200T} . These authors analyzed magnetic field data for two quasiperpendicu- 
lar shocks with moderate and high Alfven Mach number. While for moderate Ma the shock 
profiles measured by different spacecraft were approximately the same, with the exception 
of a small-amplitude wave activity in the foot, for high Ma the amplitude of the fluctuations 
attains 10 nT, making profiles considerably different for different spacecraft. However, the 
authors argued that these fluctuations stop before the ramp and do not appear to disrupt the 
shock structure; on the other hand, they didn't reject an opportunity that the fluctuations ob- 
served may represent the signatures of the unsteady shock reformation. Hereafter we report 
the first direct observation that clearly evidence the shock front reformation observed on- 
board Cluster mission on 24th of January 2001. This material was first published in |Lobziri| 
[eTST1 < [2007] >. 

The sixth problem we discuss in this Review is important for the definition of the relative 
role of dissipative and dispersive effects, namely the problem of anomalous resistivity. The 
problem of electron heating mentioned above is considered for many years to be 'solved' 
for subcritical shocks and conventional solution proposed and widely accepted is formulated 
in terms of magic words 'anomalous resistivity'. This notion introduced first by Sa gdeev| 
(1965b) and then analyzed in more detail by |Galeev| l |1976| ), who made estimates of the 
characteristic scale of the shock transition relying on ion-sound instability. Papado poulos| 
< |1985a[ > has noticed that in case of quasiperpendicular shocks the most important instabili- 
ties should be related to lower hybrid waves and has revised the model taking these effects 
into account. However there were no any measurements that might be used to confirm or 
reject theoretical models of the dissipation due to anomalous resistivity. It is worth noting 
that this problem is very important for the determination of energy redistribution between 
electrons and ions, especially for the electron heating and electron acceleration. We can not 
present here theoretical studies of anomalous resistivity, interested reader can find general 
ideas in the review papers by Galeev and Su dan] |T989| and[Papadopoulos ( 1985a). The sec- 
ond important problem where the short scale length waves are important is an energization 
of electrons within a collisionless shock. It requires the transfer of a portion of the energy 
associated with the incoming upstream plasma flow to the electron population. In order for 



14 



Krasnoselskikh et al. 



this energy transfer to occur, there has to be some media that can channel energy from the in- 
coming ion population to the electrons. One mechanism that has commonly been proposed, 
both for solar systems and particularly for astrophysical applications is based on excitation 
of lower-hybrid waves ( Papadopo uios|| 1 98 1 1?[ [Laming 2001] |Krasnoselskikh et al |1985| >. 
The increased level of electric field fluctuations in the frequency range 10 2 — 10 Hz ob- 
served in the vicinity of a quasiperpendicular shock front is usually attributed to either ion- 
sound or whistler waves. One of the most comprehensive studies of the plasma waves in this 
frequency range was conducted by |Gurnett| l |1985| l. It's main conclusion was that waves ob- 
served above local electron cyclotron frequency are Doppler shifted ion-sound waves whilst 
those below are whistler mode waves. 

The main reason of the lack of measurements of ion sound and lower hybrid waves is re- 
lated to technical difficulties of small scale electric field measurements. Recently two papers 
were published that represent first attempts to create the experimental basis for the anoma- 
lous resistivity studies in collisionless shocks. We present in this Section a short summary 
of the results obtained following to |Balikhin eFaT] ( |20031 > ; | Walker et al.| ( [2008] >. The rapid 
changes that are observed in the plasma at the front of a supercritical, quasi-perpendicular 
shock and described in previous Sectons lead to the creation of multiple free energy sources 
for various plasma instabilities. Twin satellite missions, such as ISEE or AMPTE, have pro- 
vided data for a number of comprehensive surveys of the waves observed in the frequency 
range (10~ 2 — lO^^Hz) of the plasma turbulence encountered at the shock front. The use of 
multisatellite data for wave identification and turbulence studies is limited to the analysis 
of those waves whose coherence lengths are of the same order of magnitude as the satel- 
lite separation distance. Plasma wave modes such as the ion-sound or Lower-Hybrid, that 
are supposed to play an important role at the shock front, possess coherence lengths that 
are very short in comparison with any realistic satellite separation distance l |Smirnov and| 
Vaisberg 1987). For majority of these waves the coherence length is either comparable to 
or a few times greater than their wavelength. In such cases the waves observed by differ- 
ent satellites in a multisatellite mission will be uncorrelated. This will make it impossible 
to apply wave identification methods based on intersatellite phase delays (Bali khin et al.| 
|1997a|[2003| > or k-filtering (Puujon and Glassmeier 2008 1. Nevertheless the identification of 
waves with short wavelengths and study of their dynamics remains very important because 
of their potential role in the transfer of energy associated with the upstream directed motion 
into other degrees of freedom. In the classical model of a quasiperpendicular low /3 shock 
anomalous resistivity occurs due to ion-sound turbulence in the shock front (Galeev 1976 1. 
Lower hybrid waves also may play an important role at the shock front since they also can 
be involved in resonance interactions both with electrons and ions and so may be extremely 
efficient at channelling the energy exchange between the two spices. In order to assess the 
importance of ion-sound, lower hybrid and other waves with relatively short wavelengths 
within the plasma dynamics of the shock front the mode of the observed waves should first 
be identified. The strong Doppler shift that results from the large values of wavevector |fc| 
precludes the reliable use of the observed frequency for correct identification as was done in 
many previous studies. Here we show that the data from a single spacecraft can be used to 
determine propagation modes of waves observed in the frequency range 10 2 - 10 4 Hz at the 
front of the terrestrial bow shock. A similar approach has been used by |Tjulin et al.| f2003l 
in a study of lower-hybrid waves in the inner magnetosphere. 

The lower-hybrid wave is an electrostatic plasma wave mode whose plasma frame fre- 
quency is in the vicinity of the lower- hybrid resonance frequency oj//, ~ \J Q ci Q ce where 
Q c i and f2 ce are the ion and electron gyrofrequencies respectively. The wave has linear po- 
larisation and propagates almost perpendicular to the magnetic field (cos(^s) ~ \Jm e /mi ~ 
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89°). The maximum growth rate "/max occurs when ku/k ~ LUpi/ujpe, where ku is parallel 
component of the wavevector. Since the waves are propagating in a plasma that is moving 
with respect to the satellite, their frequencies will be Doppler shifted in the spacecraft frame. 
The magnitude of this shift can be estimated using the resonance condition of the Modified 
Two Stream Instability (MTSI) IVaMaJk, = ui^. This gives a maximum estimate for the 
correction in observed wave frequency due to the Doppler shift kV sw ~ w;/j/2, here V aw is 
the velocity of the solar wind supposed to be equal to normal component of the upstream 
flow velocity. 

Current models of wave turbulence that determines anomalous resistivity in the front of 
quasiperpendicular shocks are based on the occurrence of lower hybrid waves at a shock 
front being generated due to counter-streaming populations of ions and relative motion of 
reflected ions ( |Leroy et al.|1982] > and bulk electrons and ions at the front via the modula- 
tional two-stream instability (MTSI) (Papadopoulos 1985b; Matsukiyo and Scholer 2006]l or 
modified Buneman instability. These models are often used to explain the electron acceler- 
ation observed at various astrophysical shocks such as supernova remnants (Laming 2001| >. 
However, there is currently no substantial experimental evidence that these waves do in- 
deed exist in the fronts of supercritical, quasiperpendicular, collisionless shocks. The results 
of data analysis from the Intershock electric field experiment, in which wave activity was 
observed at frequencies of a few Hertz, has been used to argue for the existence of lower hy- 
brid waves ( Vaisber g et al.|1983| . An alternative explanation, however, has been proposed in 
which Intershock may have simply observed the electric field component of whistler wave 
packets propagating in the foot region ( Krasnoselskikh et al.|199l||Balikhin et al.| [T997a; 
Walker et al. 1999a). Electric field observations of Comet Halley also showed evidence for 
waves observed in the vicinity of the lower hybrid frequency (Kli mov et al.|19 86 ). However, 
their exact wave mode was not determined. The natural way to differentiate these modes is 
to examine their polarisations. Whistler mode waves are elliptically polarised whilst lower 
hybrid waves, as mentioned above, are linearly polarised. We present here the summary of 
direct observations of the ion-sound (Balikhin et al. 2005) and lower-hybrid waves (Walker 
et al. 2008| that we complet by estimates of characteristic effective collision frequencies. 
The paper is organized as follows. In Section|2]we present the statistical studies of quasiper- 
pendicular shock ramp widths. Section [3] is dedicated to electric field scales of the ramp 
of quasiperpendicular shocks. In SectionHlthe results of evaluation of the Poynting flux of 
oblique whistler waves upstream of the shock front are outlined. Recent results on electron 
heating scale at High Mach number quasiperpendicular shocks are presented in Section[5] In 
Section[6]we use the data of measurements of lower hybrid and ion sound waves intensities 
to evaluate the characteristics of anomalous resistivity aiming to determine its role in the 
shock front formation. In Section[7]we present results of direct observations onboard Clus- 
ter satellites of nonstationarity and reformation of high-Mach number quasiperpendicular 
shock. In Section [73] we resume the results of experimental observations and discuss the 
conclusions that follow from Cluster and THEMIS observations. In Appendix[A]we present 
some short comments concerning comparison of computer simulations with the observa- 
tions. Appendix [E] contains a notation table defining the parameters used in this paper. 

2 Statistical studies of quasi-perpendicular shocks ramp widths 

As was discussed in the Introduction the characteristics of the major transition within the 
shock in which the flow deceleration and the magnetic field and electrostatic potential vari- 
ations take place are determined by the interplay between nonlinearity, dispersion and dissi- 
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Fig. 2 The total magnetic field strength as a function of distance through five low beta, quasi-perpendicular 
shocks in order of increasing ratio of criticality. The Mach number, /3, and 9g]\f for each shock are displayed. 
The data are shown at the highest temporal resolution available. The sampling rate for the first three shocks 
is 16 Hz and 4Hz for the last two (Adapted from Farri s et al.|1993) . 



pation. The presence of a population of reflected ions makes it difficult to construct a reliable 
theoretical model based on an analytical or semi-analytical description. However, the estab- 
lishment of the scales of this transition allows one to determine the characteristics of the 
dominant physical processes in play. Single satellite missions provide very poor possibil- 
ities for the reliable identification of the shock width and evaluation of the characteristic 
scales of structures within it. In such cases the spatial size of the foot or overshoot have 
been used (Balikhin et al. 1995 1 to evaluate the thickness of the ramp region. The first shock 
crossings by two satellites were studied in the frame of ISEE and AMPTE projects. These 
missions provided the first insight into the thickness of the shock transition. The decrease 
of the shock width and consequent increase of gradients as a function of increasing Mach 
number was clearly demonstrated by |Farris et"aL"| ( |1993| >. In their study the Mach number 
was normalized to the critical Mach number that determines the transition from sub-critical 
to supercritical shocks. Figure [2] (From Farr is et al.|1993[ l shows magnetic field magnitude 
measurements made by ISEE 1 for five low beta quasiperpendicular shocks ordered by the 
ratio of their Mach number to critical Mach number increasing from top to bottom. The 
increase in the gradient of the shock transition is clearly independent of the differences in 
the angles between the normal to the shock front and the magnetic field and of the value of 
beta. 

Several studies have been dedicated to the investigation of structural elements of the 
shock front making use of ISEE 1,2 magnetometer data. Scudder and co-authors (Scudder 
jetaT 1986a b c) carried out the detailed study of the shock crossing on 7th November 1977. 
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This is presumably the most detailed study of a single event, in which all the elements of the 
structure were put together and compared with detailed measurements of the particle distri- 
bution functions. These authors succeeded in relating the evolution of the ion distribution 
function to the characteristic features of the magnetic field structure and in the determina- 
tion of the major macro-features of structure of the shock front. This study concluded that 
the size of the magnetic field transition was determined by the dissipative process related to 
reflection of ions. Twin satellite measurements by ISEE provided the first indications that 
some shocks have quite narrow fronts (Newbury and Russell 1996; Newbury et al. 19981. 
Newbury and Russell ( 1996 ) reported the observation of an extremely thin, quasiperpendic- 
ular shock whose ramp width was determined to be 0.051/^ (where L t is ion inertial length), 
i.e. of the order of electron inertial length. Figure [3] shows this particular shock crossing 
(bottom panel) together with a second shock, observed under similar solar wind conditions, 
whose ramp width was 0.89Li. 

Cluster and THEMIS have provided new opportunities for a comprehensive study of 
the shock ramp scales. Recently, there have been two papers dedicated to studies of the 
ramp scales of the magnetic field transition (Hobara et al. 2010; Mazelle et al. 20101. They 
have used slightly different definition of the ramp thickness and scale and applied different 
methodology, however they came to very similar conclusions about statistical characteristics 
of the scales of the ramp transition. Hereafter we present the summary of results reported in 
these papers. 

The article by Hoba ra et al.|p010| is devoted to statistical studies of the magnetic field 
spatial scales in the ramp region of the shock front based on Cluster and THEMIS observa- 
tions. Due to their highly inclined orbit, the Cluster satellites enable the observation of shock 
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crossings away from ecliptic plane. These shocks typically exhibit Mach numbers that are 
in the lower range of the whole space of terrestrial Mach numbers, since the shock nor- 
mal deviates from the sunward direction. To increase the range of available Mach numbers, 
THEMIS shock crossings were added to the set of Cluster observations. Magnetic ramps 
cannot be always treated as uniform. Nonlinear substructures have been observed and re- 
ported within the ramp in several cases (e.g. Bali khin et al.|2 002; Wal ker et al.|20 04). The 
study of spatial-temporal characteristics of such substructures requires at least two point 
measurements separated by a distance that is sufficiently smaller than the inter satellite dis- 
tances of both THEMIS and Cluster missions, thus the authors restricted themselves with 
the study of the ramp spatial scales only. 



2.1 Criteria for Choosing Shocks and Definition of Notions "Size" and "Scale" 

Hob ara et al.| ( |2010| l have used the data from Cluster and THEMIS for a statistical study 
of the spatial size of the ramp. Both these missions assembled a huge stockpile of shock 
crossings. These data sets complement each other because of the difference in the orbits of 
Cluster and THEMIS satellites. The THEMIS orbit is close to equatorial plane providing an 
opportunity to sample the terrestrial bow shock in the vicinity of the subsolar point. Cluster 
crossings of the terrestrial bow shock occur mainly on the flanks. The solar wind flow in the 
vicinity of the terrestrial orbit is almost along Sun-Earth line, so that the Mach numbers of 
flank shocks are relatively low due to the greater deviation of the local shock normal from the 
sunward direction. Therefore, the combination of THEMIS and Cluster crossings allowed 
to cover a greater dynamical range of Mach numbers available for the analysis than each 
of these missions provides separately. Cluster crossings for two time intervals, February- 
April 2001 and February-March 2002, were used in the study. THEMIS shock crossings 
included in the study took place from the beginning of July 2007 to the end of August 2007. 
The magnetic field data used in the present paper came from Cluster and THEMIS fluxgate 
magnetometers (FGM) (Balogh et al.|1997||Auster et al.|2008) . Another reason to use the 
THEMIS data from the initial phase of the mission was that THEMIS C and D spacecraft 
separation was not very large. 

The set of shocks that have been used by Hob ara et al. 1(2010) for the study of statistical 
properties of the ramp width and gradient scale in the paper included 77 individual crossings 
of the terrestrial bow shock (30 by THEMIS and 47 by Cluster). In order to determine the 
spatial scale of the shock ramp by means of transformation from temporal to spatial vari- 
ables an estimate of the relative shock spacecraft velocity along the shock front normal was 
used. These estimates are very sensitive to shock normal definition. Thats why to perform 
reliable identification of the local normal to the shock front one of them ( Hobara et al . 20 1 1 
compared normals found making use of four different methods, using the model shape of 
the terrestrial bow shock similar to |Farris et al.| ( |1991[ l, using timing differences methods 
between the 4 Cluster spacecraft shock crossings, minimum variance and coplanarity the- 
orem. In order to validate the results the evolution of the magnetic field component along 
the normal direction B„ was used. The reliability of the shock normal identification served 
as the only shock selection criteria. Those shock crossings, for which the calculated normal 
could not be considered reliable (because of the B n evolution or due to large discrepancy 
in the shock normal directions found by different methods), have been excluded from con- 
sideration. The relative shock spacecraft velocity V ss has been calculated using the shock 
normal direction, satellite separation vectors and time delay between two subsequent shock 
crossings. 
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The second recent study dedicated to the same problem ( |Mazelle et al. 2010) was also 
based on Cluster magnetic field measurements during spring seasons of 2001 and 2002 cor- 
responding to small interspacecraft separation (100 to 600 km typically). The shock param- 
eters (angle between upstream magnetic field and local shock normal 8g n , Alfvenic Mach 
number A/4 and ion beta /%) were computed from the data of Cluster and ACE. 

Mazelle et al. (2010) selected the shocks for the statistical study according to follow- 
ing criteria. First, in order to restrict the study by almost perpendicular shocks the shocks 
with Obu as close as possible to 90° were chosen. The other criteria were the existence of 
well-defined upstream and downstream intervals for the 4 s/c, the stability of the upstream 
averaged field from one s/c to another, the validity of the normal determination by checking 
the weak variability of the normal field component B n around the ramp and low value of B n 
upstream for 6g n to be close to 90°. Only 24 from 455 crossings of Cluster satelite quartet 
in 2001 and 2002 were left with all criteria validated. This means that 96 individual shock 
crossings were analyzed. Selected 6g n values were chosen to be in the range from nearly 
90° to 75° but about 80% of the shocks selected were above 84° . Mach numbers Ma were 
found to be equally distributed from 2 to 6.5 and corresponding between very low values 
to 0.6 but with 67 of values less than 0. 1 . 

The major difference between two works consists in using different methods of the de- 
termination of the shock ramp thickness. If the beginning of the ramp region was defined 
quite similarly as the beginning of the monotonous increase of the magnetic field, the end of 
ramp or exit from the ramp region was determined differently. |Hobara et al.| ( [2"010| ) defined 
the ramp crossing duration as a time interval between the upstream edge of the ramp and the 
maximum of overshoot. The spatial size (width) of the magnetic ramp has been estimated 
as a product of V ss and the duration Dt of the ramp crossing. Mazel le et al.| ( |2010[ l deter- 
mined at first a stationary asymptotic level of the magnetic field that might be considered 
as satisfying to Rankine-Hugoniot conditions. To this end a downstream interval where the 
magnetic field magnitude is quite steady was used. It is then considered as giving an approx- 
imate estimate of the value of the magnetic field corresponding to exit from the ramp/entry 
in the overshoot. From the initial values of the entry in and exit from the ramp a linear fit 
of the data points inside the estimate time interval is made and this later one is allowed to 
vary. The choice of a linear fit allows excluding any pollution of the ramp region by a part 
of the extended foot. The same analysis was repeated for all four satellites for each shock 
crossing. The steeper slope found for the ramp defines the 'reference satellite'. The times of 
the middle of the four samples of ramps for one shock crossing are then used to compute 
both the shock normal and velocity in the GSE frame by the 'timing-method' described in 
|Horbury et al.| (2~002[ ). This makes it possible to derive the 'apparent' width (along each 
satellite trajectory) and compare between the 4 s/c. Then, the velocity vector of the shock 
in each s/c frame is computed. Its angle with the shock normal allows reconstructing a local 
profile along the normal and determining the real local spatial width of each shock sub- 
structure. Readers interested in more technical details can find them in Maz elle et al.| ([2010l 
and |Hobara eTaT]p010l l. 

A scatterplot of the spatial sizes of the ramp as a function of Alfven Mach number 
is shown in the left hand panel of Figure [4] The lefthand panel shows the scale sizes in 
terms of the ion inertial length whilst the right hand panel shows the width in terms of 
the electron inertial length. The width of the magnetic ramp varies by about an order of 
magnitude from L r = 1.4L^ (« 60L e ) to O.lLi (w 4L e ). The general trend in Figure]?] 
indicates that the magnetic ramp becomes thinner with increasing of Alfvenic Mach number. 
This trend is evident even without taking into account two shock crossings with peculiarly 
high Mach numbers M in the range 17-20 that correspond to the two markers in the bottom 
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Fig. 4 Scatterplot of experimentally derived shock size (left panel) and shock spatial scale (right panel) 
normalized to the ion inertial scale length c/oj p ; (left axis) and electron inertial scale c/ui pe (right axis) as a 
function of Alfven Mach number. The dashed line represents the averaged values of shock width averaged 
over shocks with Alfven Mach number in the ranges 2-4, 4-6, 6-8, 8-10, 10-12. The vertical lines represent 
the statistical error bars for each of these Mach number ranges. (Adapted from Hoba raet al.|2010) 



right corner of the scatterplot. The figure also shows a distinct decrease in the maximum 
width of the shock front with increasing Mach number. To make these tendencies more 
clear the characteristic width of the magnetic ramp averaged for the shocks with Alfvenic 
Mach numbers in ranges 2-4, 4-6, 6-8, 8-10 and 10-12 (dashed curve) are presented. 

The vertical lines on Figure [4] (left panel) represent the statistical error bars for each 
range of Mach numbers 2-4, 4-6, 6-8, 8-10 and 10-12. The decrease of statistical errors 
with the Mach number is in complete agreement with the significant decrease of the maxi- 
mum shock width, while the minimum shock width undergoes much smaller changes. 

The right hand panel of Figure [4] shows the scatterplot of the spatial gradient scale of the 
magnetic ramp. It clearly demonstrates the same characteristic features as were evidenced 
in the left hand panel for the ramp width, namely a quite wide range of values, especially 
for low Mach number shocks and the trend toward shorter scales with the increase of the 
Mach number as well as the decrease of the maximum gradient scale while Mach number 
increases. As the change of the magnetic field for all chosen shocks exceeds upstream mag- 
netic field Bo (for many of them quite significantly) the values of the gradient scale are 
smaller then the width of the corresponding shocks. The ramp gradient spatial scale varies 
in the range 0.05 - 0.82L i (2 - 35L e ). 



2.2 Statistical Analysis ( |Mazelle etal.|2010[ l 

Figure[5]displays the thinnest ramp found among each quartet of crossings for each individ- 
ual shock versus 8g n . There were no simple relation found. However, it is unambiguously 
established that many observed thinnest ramps are less than 5 c/cj pe thick and there was 
found an apparent trend for lower values as 8g n — > 90° . 



Dynamic quasiperpendicular shock 



21 




Bn 



Fig. 5 Thinnest ramps observed versus shock 
9g n . Reprinted with permission from (Mazelle 
|et al.|2010> . Copyright 2010 , American Institute 
of Physics. 
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Fig. 6 Histogram of the 96 shock ramp thick- 
nesses. Reprinted with permission from I Mazelle 
|et al.|2010) . Copyright 2010 , American Institute 
of Physics. 



The histogram of all ramp thicknesses in Figure [6] reveals the predominance of narrow 
ramp width with a Gaussian-like regular decrease towards an asymptotic limit that is still 
less than c/uj p i. 

The authors came to conclusion that their analysis confirms statistically that the mag- 
netic field ramp of the supercritical quasi-perpendicular shock often reaches a few c/uj pe . 

So, the results of two independent studies by two different groups come to the same 
conclusion, the ramp width for quasiperpendicular high Mach number shocks as seen in 
magnetic field is of the order of several c/u! pe and is in a perfect agreement with estimates 
determined from the dispersive scale of a nonlinear whistler wave, modified by the presence 
of reflected ions l |Galeev et al.|1988|[1989l|Krasnoselskikh et al.|2002) . 



3 Electric field scales of the ramp of quasiperpendicular shocks 

As mentioned in the introduction, there have only been a few reports regarding the scales and 
structure of the electric field transition at quasi-perpendicular shocks. Based on laboratory 
experiments, in which the conditions are not exactly the same as in space plasmas, Ese levich| 
|et al.| ( [797T] l reported that the major change in potential across the shock occurs within the 
magnetic ramp region. 

Figure[7]is a sketch (based on the results of Eselevich et aL] ( [T97T] l) of the change in the 
magnitude of the magnetic field and the accompanying change in the electrostatic potential. 
These authors interpreted it as a viscous subshock similar to isomagnetic jump. From Fig- 
ure[7]it can be seen that there are two different length scales that may be associated with the 
change in the electrostatic potential as the shock is crossed. The first, indicated by the lightly 
shaded bar at the foot of the plot, shows that overall the potential changes on scales similar 
to that of the magnetic ramp region in agreement with the results of Lembe ge et al.| fl999). 
This corresponds to an enhancement of the electric field observed as the shock is crossed. 
The second scale, indicated by the darkly shaded bar, corresponds to a region within the 
shock front in which a large increase in the potential is observed over a small spatial scale. 
Such changes in the potential result from large amplitude spike like features in the electric 
field. 
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Fig. 7 Sketch of the changes observed in the magnetic field and electrostatic potential during the crossing of 
a quasi-perpendicular shock (based upon the experimental results of I Eselevich et al.|197T) ) 



The results reported here present a study of the large amplitude, short duration features 
in the electric field observed by the Cluster satellites during a number of crossings of the 
quasi-perpendicular bow shock published in Walker et al. (2004). These features contribute 
significantly to the overall change in potential observed at a shock crossing but their short 
duration implies that they are very localised. The aim of the study was to determine their 
scale size and amplitudes. These parameters were studied in relation to the upstream shock 
parameters. 

A total of 54 shock crossings, occurring on 1 1 separate days were investigated but not all 
could be analysed fully for various reasons such as unavailability of certain data sets, or the 
accuracy of the shock normal determination. In this section we present a case study of the 
electric field within the shock front, namely the crossing that occurred on March 31 st , 2001 
at around 1718 UT. On this day the conditions in the solar wind were some what abnormal 
due to the passage of a CME. Measurements in the solar wind by Cluster indicated that the 
magnitude of the magnetic field was of the order of 30nT, the normal for this shock (based 
upon FGM crossing times) is ng = (0.94, —0.17, 0.293) (in the GSE frame), and the shock 
velocity was determined to be 48.92kms~ 1 . Other relevant parameters are m 87°- 
and a density n ~ 19cm~ 3 . The high value of the field resulted in an unusually small 
P ~ 0.07. The Alfven Mach number for this shock (A/4 ~ 3.6) lies close to the First Critical 
Mach number and to the whistler critical Mach number so the conditions of the solar wind 
are quite favourable for the formation of quasi-electrostatic sub-shocks at the shock front 
l |BaUkhiii et al.[2002) . Figure [gj (adopted from Wal ker et al.] ( |2004| ) shows an overview of 
the magnetic and electric field measurements made by FGM and EFW respectively during 
this shock crossing. The top panel shows the magnitude of the magnetic field measured by 
FGM. Initially, all four Cluster spacecraft are in the solar wind just upstream of the outward 
moving bow shock which subsequently swept over the satellites in the order C4 (17: 17:43.5), 
C2 (17:17:45.5), CI (17:17:48.5), and finally C3 (17:17:53.5). The magnetic field profiles 
show a set of clean shock crossings that possess clearly discernible foot, ramp and overshoot 
regions. The second panel shows the magnitude of the electric field measured by EFW in 
the spin plane of each satellite (\E\ 2 = E% + Ey). In the solar wind, the typical magnitude 
of the electric field is around 14mVm _1 in the satellite spin plane. It is possible to estimate 
the E z component of the upstream electric field assuming that E • B = 0. This assumption 
is valid for estimating the field upstream and downstream of the shock but not within the 
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Fig. 8 Overview of the shock crossing on March 31 st , 2001 at 1718 UT. The top panel shows the magnitude 
of the magnetic field measured by FGM. The second panel shows the magnitude of the electric field measured 
in the satellites spin plane. The lower two panels show the spin plane components E x and E y . (Adapted from 
|Walker et al."|2004t 

shock region itself. Upstream of the shock, Ez » 5mVm _1 . This value is higher than the 
measured E x component (~ 2.5mVm _1 ) and less than the E y component (-13mVm _1 ). 
Comparing the top two panels it can be seen that the disturbances measured in the electric 
field begin in the foot region of the shock and continue until the satellites are downstream 
of the overshoot/undershoot. These general disturbances have amplitudes generally in the 
range 5 — 30mVm _1 . During their crossings, each of the satellites recorded a number of 
large amplitude, short duration features in the electric field. The largest of these spikes have 
maximum amplitudes of approximately 30, 40, 60, and 65mVm _1 for satellites 1, 2, 3, and 
4 respectively above the field measured in the solar wind just upstream of the shock front. 
These values represent lower limits of the strength of the electric field since the component 
perpendicular to the spin plane is not considered. They are seen to occur within the ramp 
region but there is no strong feature within the FGM data with which they correlate. It is 
also observed that the largest electric field peaks observed on each satellite appear to occur 
in pairs which may suggest field rotation. The two lower panels show the components of the 
electric field measured in the satellite spin planes. Both panels show that the components of 
the field exhibit a twin peaked structure, similar to that observed in the field magnitude and 
that the direction of rotation is the same for both peaks. Thus the overall structure is not due 
to a single rotation of the field. Our goal is a statistical study of the widths of these short 
living, large amplitude features. 

Using the four point measurements one can determine the occurrence time of these peaks 
in the electric field and hence compute a normal. Examining the Ex component, the time dif- 
ferences between the observations of the first peak in the electric field are Ati2 = —3.01s, 
Atis = 5.03s, and Atn = —5.35s. When coupled with the respective positions of the 
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Fig. 9 The FGM magnetic and EFW electric fields measured by Cluster 1 on March 31 st , 2001 around 17:18 
UT. The magnetic field magnitude is shown by the magenta line. The spin plane electric field magnitude, and 
E x and E y components are shown in red, blue and cyan respectively. The yellow regions highlight the periods 
when large amplitude short duration spikes in the electric field are observed. The black line (Y scale of RHS) 
represents the change in potential within the shock. (Adapted from Walke r~et al.|2004) 

satellites this yields a normal direction rig = (0.946,-0.155,0.283) and a velocity of 
~ 50kms~ 1 . The difference between this normal n e and that determined from the magnetic 
field (rifl) is less than a degree. Thus it appears that the electric field spikes correspond to 
layers within the overall shock structure. Figure|9]shows the results from Cluster 1 in greater 
detail. The magenta line shows the magnitude of the magnetic field. The foot region was en- 
tered around 17: 17:47.3 UT whilst the ramp was crossed between 17:17:48.3 and 17:17:48.9 
UT. Several large spikes in the electric field are observed in the region of the foot and shock 
ramp. The three largest occur around 17:17:48.2 (20mVm _1 ), 17:17:48.5 (30mVm _1 ), and 
17:17:48.6 (15mVm _1 ). Their short duration implies that their scale size is of the order 
3 — 5c/uip e . The black line in Figure |9]represents an estimation of the electrostatic potential 
measured in the normal direction. This was calculated by removing an average of the field 
measured in the region just upstream of the shock from the field measured within the shock 
region and then integrating the projection of this electric field along the normal direction. 
Whilst the actual potential cannot be calculated due to the incomplete vector measurements, 
it can be estimated by assuming that the field perpendicular to the spin plane E z = 0. This 
assumption is valid because for this particular shock, the normal lies very close to the spin 
plane. This calculation can be used to show that the largest jumps in the potential coincide 
with the spikes observed in the electric field and that these occurrences contribute a signif- 
icant fraction of the total potential change observed at the shock. During this period, the 
electric field enhancements contribute around 40% of the total change. 




3.1 Scale size 

The preceding sections have presented evidence for localized increases in the electric field 
strength measured as the satellite traverses a quasi-perpendicular bow shock. All shocks an- 
alyzed show evidence for an enhancement in the background electric field. In most cases, 
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Fig. 10 Histogram of the scale sizes for the spike-like enhancements observed during a number of crossings 
of the quasi-perpendicular bow shock. (Adapted from Walker et al. 2004) 

the region in which this field enhancement occurs lasts longer than the crossing of the mag- 
netic ramp. The field typically increases of the order l-3mVm -1 above that measured in the 
solar wind. However, as has been noted above, the turbulence in this region is dominated 
by spike-like fluctuations lasting a few milliseconds and with magnitudes of typically 4- 
20mVm~ 1 with a maximum magnitude of the order of 70mVm _1 . This existence of large 
gradients in the electric field has repercussions for processes involved in the heating of elec- 
trons. In the presence of strong electric field gradients the electron gyration frequency can 
deviate from its classically calculated value (Cole 1976; Balikhin et al. 1998), leading to 
an increase in its Larmor radius and the possibility of a breakdown in adiabaticity l |Cole| 
|1976 ; Balikhin et al.|1998) . Having shown that the spikes observed in the electric field at 
the front of a quasiperpendicular shock appear to be physical structures that form a layer 
within the shock front as opposed to being the result of noise in the data or motion of the 
shock a statistical study of these features was performed to investigate their relationship to 
the properties of the shock front. Now we shall present statistical study of the data collected 
from a number of such spike-like features. 

Figure [TO] shows the distribution of the scale sizes determined from the event duration 
and the shock velocity evaluation of these features in terms of the electron inertial length. 
The scale size of these events will be unaffected by the incomplete vector measurements of 
the electric field. The vast majority of crossings have scale sizes of the order of 1 — 5c/oj pe - 
The data that form tail of the distribution at longer scale sizes typically comprise events 
that have a multi-peak structure. This type of event represents an upper limit to the scale 
size of these short-lived events. In comparison, the typical scale of the magnetic ramp is 
characterized by the ion inertial length (Newbury and Russell 1996) although these authors 
also report one particular shock as having a ramp scale as small as 0.05c/uj pi or 2c/w pe - 

Figure [TT| shows the relationship between the Mach number and scale size of the spikes 
observed in the electric field. From the figure, it is clear that the scale size has a lower 
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Fig. 11 Dependence of scale size on upstream Mach number. (From Walker et al. 2004 



limit that decreases as the Mach number increases. One should notice that these results 
represent the tendency rather than the proof, the number of points at high Mach numbers is 
not sufficient for valuable statistical study. 

Figure [12] shows a scatter plot of the relationship between 8g n and the scale size of the 
electric field enhancements. In general there appears to be a broad range of scales. However, 
as 8g n approaches 90°- the scale length decreases. For the shocks analysed with 8g n close 
to 90° - the scale lengths are of the order of 2c/cj pe . This compares favourably with theoret- 
ical estimates that for shocks close to perpendicular the scale width is estimated to be of the 
order of the electron inertial length as proposed by Karpman|( |1964| l; |Galeev et al.| fl988). 
This tendency corresponds exactly to dispersion dependence upon the angle 8^ n . 



3.1.1 Amplitude 

The examples presented above show that the increase in the electric field (AE = E spi ], e — 
Eupstream) observed during encounters with these spike-like structures varies between 4 
and 70m Vm -1 above the average field that is measured in the solar wind just upstream of 
the shock. In this section the relationship between this change (AE) and the shock Mach 
number, and the angle 8g n is presented. Figure[T3]shows a scatter plot of the peak amplitude 
observed in the electric field spike event (AE) as a function of the shock Mach number Ma- 
For shocks whose Mach number Ma > 5 there is a fairly constant trend in which AE < 
15mVm _1 . In the Mach number range 3 < Ma < 5 the range of observed amplitudes 
varies between 5 and 60mVm~ 1 . It appears that in this Mach number range the electric 
field potential becomes more important than for low and high Mach number shocks that 
corresponds to dependence of the electrostatic potential upon the Mach number ( Di mmock| 
et al.||2012| . The red crosses highlight the shocks observed on March 31 st , 2001. All of 
these shocks fall into this range of Mach numbers. This set of shocks seems to possess 
Mach numbers corresponding to supercritical range and having large number density, at 
about 7%, of alpha particles (Maksimovic, private communication). Their structure seems 
to resemble that of electrostatic sub-shocks similar to those observed in laboratory plasmas 
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Fig. 12 Dependence of scale size on 9g n . (From Walker et al. 2004 



(Eselevic h et al.|1971| l. A characteristic signature of sub-shocks is the occurrence of small 
scale electrostatic fluctuations such as those observed on this particular day. Ion sound sub- 
shocks have been observed in laboratory plasmas with scales of the order of 100 Debye 
lengths. For the shocks observed in March 31 s * 2001, the scale is closer to characteristic 
scale of the fast magnetosonic mode (Balikhin et al. 2002]l. 

The relationship between AE and is shown in Figure [14] It clearly shows that 
as Q^n approaches 90°- the range of the observed amplitudes of the electric field spikes 
increases. 



3.2 Conclusions 



In this Section we presented the changes observed in the electric field during the crossing of 
a number of quasi-perpendicular bow shock. It is shown that the electric field is enhanced 
during the crossing of the shock and that the scale size over which this enhancement is ob- 
served is larger than that of the macroscopic magnetic ramp region. Within the whole shock 
region, short lived electrostatic structures are observed that are intensified in the ramp re- 
gion. The scale size of these structures is of the order of a few c/uj pe : and was shown to 
decrease as 8g n approaches 90°- which corresponds to the dependencies following from 
theoretical model based on consideration of the shock as mainly dispersive nonlinear struc- 
ture (Galeev et al. 1988; Krasnoselskikh et al. 2002l. The amplitudes of these structures is 
typically of the order of 5-20mVm -i : but under special circumstances may reach as high as 
VOmVm" 1 . The highest amplitudes appear to be observed for shocks whose Mach number 
is in the range 3 to 5. This may be an indication that such shocks have quasi-electrostatic 
sub-shocks inside the main ramp transition. It was also demonstrated that these small scale 
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Fig. 13 Scatter plot showing the relationship between the amplitude of the electric field spikes as a function 
of Mach number. The red crosses are used to highlight the data for the shocks that occurred on March 31 st , 
2001. (Adapted from |Walker et al.|2004^ 
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Fig. 14 The relationship between AE and 0B n - (Adapted from Walker et al. 2004 



structures make a substantial contribution to the overall change in potential observed across 
the shock and that the potential change is not linear. 



4 Dispersive nature of High Mach number shocks: Poynting flux of oblique whistler 

waves 



It is well known that a subcritical shock has a nonlinear whistler wave train upstream of its 
front (Sagdeev 1966; Mellott 1985). The major transition of such a dispersive shock, the 
ramp, behaves as the largest peak of the whistler precursor wave package ( K arpman et al.| 
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[T973) |Kennel et aL][T985l |Galeev et al"][l989l |Krasnoselskikh et al.|[2002] i. The presence 
of whistler/fast magnetosonic precursor wave trains in supercritical shocks was experimen- 
tally established in |Balikhin eTaT1 ( |1997b| l; |Krasnoselskikh et al.| ( [T991) ; |Oka et al.| ( p006) . 
These whistler waves have rather large amplitudes and their role in energy transformation 
and redistribution between different particle populations and in the formation of the shock 
front structure is still an open question. The energy source responsible for the generation of 
these waves is the subject of active debate in shock physics (see Galeev et al. 1988, 1989, 
|Krasnoselskikh et al.|2002| |Matsukiyo and Scholer|2006[ |Comigel et al.|2011) . Often the 
precursor waves are almost phase-standing in the shock frame. However, if they are gener- 
ated by the ramp region as the dispersive precursor their group velocity can still be greater 
than zero in the shock reference frame, which would allow energy flow in the form of Poynt- 
ing flux to be emitted towards the upstream of the shock transition. On the other hand, if the 
waves are generated by instabilities related to reflected ions their energy flux will be directed 
from the upstream region towards the shock ramp. The goal of this Section is to address this 
problem, to present the direct measurement of the Poynting flux of the upstream whistler 
waves aiming to establish the direction of the Poynting flux. 

Below we establish the energy source of the waves by calculating the Poynting flux of 
the waves in the Normal Incidence Frame (NIF) of the shock, using multi-satellite Cluster 
data from crossings of the Earth's bowshock (Escoubet et a l.|1997[|Bale et al.|20 05). Two 
events with supercritical Alfvenic Mach numbers are analyzed. In both cases it is found that 
the shocks show dispersive behaviour with the Poynting flux directed in upstream direction. 

Poynting flux is not a Lorentz invariant and therefore depends on the frame of reference. 
To evaluate the value and direction of the Poynting flux with respect to the the shock we 
transform the electric field to the Normal Incidence Frame (NIF). The normal n = +x which 
also serves as the ^-coordinate direction in the NIF system is obtained by four-spacecraft 
timing, z is the direction of maximum varying magnetic field obtained from a minimum 
variance analysis, and y is the direction of the convection electric field which completes the 
right-handed system. 

The transformation from the spacecraft frame to the NIF is given by E j^jp = E sc + 
v x B. The total velocity required for this transformation is defined by v = v s ^ + vnif 
where v s /, = v^n. is the shock velocity, vpfjp = nx (v u x n) is the NIF velocity and v u 
is the solar wind velocity. 

A general shift of reference frame, coordinate transformation, and evaluation of the 
complete Poynting vector requires knowledge of the full six-dimensional electromagnetic 
field (three electric and three magnetic components). The Cluster spacecraft, however, only 
measures the two components of the elecric field in the spin-plane of the spacecraft, while 
the third component normal to the spin-plane is not measured. To reconstruct the third com- 
ponent we use the assumption that for the wave electric and magnetic fields the condition 
E ■ B = holds. While this is a true condition for the cross-shock (DC) electric field, it 
holds well for whistler wave electric fields at lower frequencies. 

We study two quasi-perpendicular high Mach number shocks encountered by the Clus- 
ter multi-spacecraft mission (Esc oubet et al.|[T997[ ). The first shock was observed around 
04:53:40 Universal Time (UT) on 20-Jan-2003, and the second around 07:07:00 UT on 24- 
Jan-2001. We use data from the EFW (electric field) jGustafsson eTaiJlWj) , FGM (DC 
magnetic field) (Balogh et al. 1997) and STAFF (wave magnetic field) instruments (Decreau 
|et al.|1997| from spacecraft 2 (for the 2003 shock) and spacecraft 3 (for the 2001 shock). 
The shock normal n is established by assuming a planar shock and using the time of crossing 
of the four spacecraft and their relative positions (Pasch mann and Daly|1998| . 
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Fig. 15 Magnetic and electric fields in the Nor- 
mal Incidence Frame (NIF) of a high Mach num- 
ber shock, a) Power spectra of the magnetic field 
(STAFF). The black line is the DC total mag- 
netic field, included to show the waves in relation 
to the shock ramp structure, b) Power spectra of 
the electric field (EFW). The data gap is due to 
instrumental interference, c) The magnetic field 
in NIF coordinates Bjsjjp. (Adapted from|Sund-| 
|kvist et al.|20T2> 



Fig. 16 Poynting flux in the Normal Incidence 
Frame (NIF) of the same shock as in Figure j!5( . 
a) Poynting flux S\ | projected on the local Bo in 
the NIF, where red corresponds to the upstream 
flux away from the shock, b) Angle s be- 
tween k and magnetic field Bo - The yellow line 
represents the average over all frequencies (right 
scale), c) Angle between Poynting flux S and 
B . (Adapted from |Sundkvist et al.|2012) 



The first shock analyzed had an upstream 8g n ~ 85° and an Alvenic Mach number 
Ma ~ 5.5. The electric and magnetic fields in the shock front region are characterized by 
waves, with stronger amplitudes closer to the ramp, see Figure (15) . The waves have fre- 
quencies f c p < /, where f cp ~ 0.1 Hz is the proton gyrofrequency, and right-handed polar- 
ization looking along the magnetic field vector and thus belong to the magnetosonic/whistler 
mode. The direction of the wave- vector k was determined by the Means method (Means 
1 1972) , which uses the imaginary part of the three-dimensional magnetic field spectral ma- 
trix. The angle O^b between the wave vector and the local ambient magnetic field is shown 
as a function of frequency in Figure (16)b. The average value (8 kB ) in the shock front region 
is ~ 10 — 50° (right-hand scale). The whistler waves are thus oblique with respect to the 
local magnetic field, as well as to the shock normal. The angle increases continously as the 
shock front is approached and O^g — > 90° at the ramp, reflecting the quasi-perpendicular 
nature of the shock. This smooth transition stresses the nature of the shock as a dispersive 
nonlinear whistler wave. 

Since Poynting flux is a second-order quantity the electric and magnetic fields in the NIF 
were wavelet transformed (Morlet width 5.36) and the cross-product Sf = l//xoE* x Bj 
formed in frequency space. The calculated Poynting flux is therefore distributed in both 
time and frequency. The projection of the Poynting flux distribution along the magnetic 
field S\ i 



'II 



Sf ■ Bq /| Bo | using an instantaneous value of Bo is plotted in Figure 1 16 la 
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04:53:40 
20-Jan-2003 



Fig. 17 Poynting flux integrated along the spacecraft trajectory. The blue line is a projection along the 
ambient magnetic field J S^dt and the green line is a projection o n the shock normal f S • ndt. The red line 
shows the scaled magnetic field B a for reference. (Adapted from Sundk vist et al.|2012) 



where the colors red (upstream) and blue (downstream) show the direction of the flux. We 
note that in the front region of the shock the Poynting flux is everywhere directed upstream 
(red), away from the shock. In the downstream area there is a mixture of blue, green and 
red, where there is more turbulence and the waves are no longer coherent. The upstream and 
slightly oblique direction of the Poynting flux is further quantified in the instantaneous angle 
8g b between the Poynting flux and the ambient magnetic field, plotted in Figure \\6)c. 



Figure shows the Poynting flux along the spacecraft trajectory, with integrated power 
over frequencies corresponding to the waves in Figure l|16|>a, 2 < / < 10 Hz. In this 
figure the slope is the important characteristic. Positive slope means Poynting flux carried 
upstream, and negative slope downstream. From the figure it is evident that the source of 
the Poynting flux is associated with the shock ramp. The data gap and associated plateau are 
due to instrumental interference. 

Another important characteristic established is that the Poynting flux direction is oblique 
with respect to the shock normal as well as the background magnetic field. This can be 
explained by analyzing how the phase velocity for whistler waves depends on this an- 
gle. The phase velocity of a wave propagating in the plane of the shock normal n and 
background magnetic field Bo, having an angle a with respect to the shock normal is 

V p h = \ -J ^ cos(#Bn — a) . Its projection on the direction of the shock normal is V p h,n = 



V ph cos a 




Bn " 

(#Bn — Q ) cos a. Its maximum value can be found to be equal to 

max(V p /j,n) — I y'^T" (1 + cos ^Bri)i mus the projected phase speed can be larger than the 
whistler critical velocity given above. The above analysis also explains the observation of 
oblique whistler wave trains found in computer simulations of purely perpendicular shocks 
l |Hellinger et al.|2007| >. So even in the case of shocks having Mach numbers larger than the 
whistler critical Mach number, whistler waves oblique with respect to the shock normal can 
remain quasi-standing. 

The second analyzed shock crossing on 24-Jan-2001 is shown in Figure This is a 
reforming high Mach number shock (Ma ~ 11) and has been analyzed in detail in |Lobzin| 
|et al.| ( p007) . Both of the shocks discussed by were analyzed using wavelet as well as Fast 
Fourier Transform (FFT) dynamic spectra techniques. We present the second shock using 
the FFT analysis, to show that the conclusions are not technique dependent. The upstream 
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STAFF/EFW, Cluster 3, 24-Jan-2001 07:04:57 




UT24-Jan-2001 

Fig. 18 Poynting flux derived from electric and magnetic fields for a high Mach number shock, a) Wave 
magnetic field and averaged Bq. b) Wave electric field, c) S\\ normalized by its standard deviation (yellow 
and red corresponds to upstream flux), d) S||. e) log 10 Su. f) Angle of S to Bq (red meaning upstream). 
(Adapted from |Sundkvist et al.|2012) 



whistler waves, Figure <|18|)a,b, again have an overall Poynting flux upstream, away from the 
shock in the normal incidence frame, evident from the red and yellow (upstream) colors of 
S|| (Figure ( 18 i, panels c through f). For this shock the ambient magnetic field was directed 
in the opposite direction, so that 180° (red) means upstream in Figure 1 18 1. 



The power flux given by the Poynting vector shows unambiguously that they carry en- 
ergy over a broad frequency range from the shock ramp towards the upstream solar wind, 
starting from the position of the shock front. This leads to conclusion that the results of 
the analysis are consistent with a theoretical model (Galeev et al. 1988; Krasnoselskikh 
|et al.|2002| > that considers the shock steepening to be balanced by the effect of dispersion 
in addition to dissipation. As the shock steepens, nonlinearities transfer energy to shorter 
wavelengths of the spectrum, and is ultimately carried away from the shock as dispersive 
whistler wave trains. This analysis demonstrates that for high Mach number shocks, dis- 
persive effects are dominant for the formation and stability of the shock front. Since the 
whistler waves are strongly damped upstream of the shock, we infer that they can play the 
role of an intermediate step in the energy re-partition problem, with the energy ultimately 
being dissipated through wave-particle interaction. 
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Fig. 19 Magnetic field data at a crossing of the Earth's bow shock by the 4 Cluster spacecraft on 9 Jan. 2005. 
Traces have been shifted by 20 nT for clarity. The dashed lines show the times of the steep ramp. (Adapted 
from |Schwartz et al.|20lT) . 



5 Electron Heating Scale at High Mach number Quasiperpendicular Shocks 



From the discussion in the previous sections the energy repartition amongst particle popula- 
tions in quasiperpendicular shocks is a multi-scale process related to the spatial and temporal 
structure of the electromagnetic fields within the shock layer. While the major features of 
the large scale ion heating are known, the electron heating and smaller scale fields remain 
poorly understood and controversial. In this section we will discuss the scale of the electron 
temperature gradient based on the possibility of obtaining unprecedented high time reso- 
lution electron distributions measured in situ by the Cluster spacecraft recently discussed 
by |Schwartz et al.] ( |2011[ l. The authors discovered that approximately half of the electron 
heating coincides with a narrow dispersive layer several electron inertial lengths (c/u) pe ) 
thick. Consequently, it gives one more argument that the nonlinear steepening is limited by 
wave dispersion. The DC electric field associated with the electron pressure gradient must 
also vary over these small scales, strongly influencing the efficiency of shocks as cosmic ray 
accelerators. 



The 4 Cluster spacecraft (Escoub et et al.| 1997] ) are unique in their ability to remove the 
time-space ambiguity in time series data taken by in situ space plasma instrumentation. By 
timing the passage of an event at each corner of the tetrahedron formed by the 4 spacecraft, 
the planar orientation and speed of the event can be determined. We employ this technique to 
convert the time series of data to distance along the shock normal Schwartz ( 1998 ). Figure[l9] 
illustrates the identification of the steep shock ramp that we use as event times. 

The electron instrument on Cluster measures fluxes at several energies in a half-plane 
containing the spacecraft spin axis. These measurements form an azimuthal wedge divided 
into 12 polar directions from aligned to anti-aligned with the spin axis, and are repeated 
at 125-250 ms intervals. A full 3D distribution covering all azimuths is thus built up over 
1 spin (~ 4 s). However, when the magnetic field is roughly aligned with the spin axis, 
each wedge contains a full set of pitch angles from 0° to 180°. Under these circumstances, 
and assuming gyrotropy, the full pitch angle distribution function is available at < 250 ms 
resolution. 
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Fig. 20 Overview of data from Cluster 2 on 2005 Jan 9. From top to bottom: Omni-directional electron 
energy-time spectrogram @ 250 ms resolution, electron pseudo-density, electron pseudo-temperatures (see 
Method), magnetic field magnitude, and field components. Arrows in the fourth panel show locations of the 
cuts presented in Figure [22] (Adapted from Schwartz et al. 201 1 1. 



We rebin the raw electron data into pitch angles a relative to the instantaneous magnetic 
field. We calculate pseudo-densities and temperatures for each pitch angle bin as if the dis- 
tribution were isotropic, e.g., n(90°) = A-k J f(v, a = 90°) v 2 dv. These pseudo-moments 
better characterise the phase space distributions in the ||,_L directions than the full Tii j_ 
moments (cf. Fig 9 of |Mitchell et al.|2012| >. 



5.1 Results and Conclusions 

An overview of the data for 2005 Jan 9 is shown in Figure|20] The transition from unshocked 
solar wind plasma to the shocked magnetosheath occurs around 22:15:30. Although the 
solar wind flow is a factor of 10 slower than the electron thermal speed, some residual 
modulation at the spin period is evident in the data. We have averaged the parallel and anti- 



parallel (a = 0, 180°) moments so that the second and third panels of Figure 20 reveal the 
pseudo-parallel and perpendicular moments. Note that the pseudo-densities n(a) are not, 
and from their definition above need not be, equal. The bottom two panels show increasing 
oscillations and a gradual "foot" ahead of a steeper magnetic "ramp" region. The dominant 
z magnetic field component is nearly aligned with the spin axis, enabling the parallel and 
perpendicular moments to be available in every 0.25 s wedge as described above. Figure [20| 
already suggests the main result namely that the rise in electron temperature closely follows 
even the steepest ramp of the magnetic field. 
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Cluster 2 2005-01-09 22:15:30 
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Fig. 21 Magnetic field (solid) and electron temperature (symbols) as a function of distance from the shock 
ramp. Roughly half the temperature rise occurs within the region 17.3 km wide between the das hed vertical 
lines corresponding to 6.4 electron inertial lengths (c/ui pe ). (Adapted from Schwartz et al.|20iT|. 



Figure [2TJ shows that both the parallel and perpendicular electron temperatures closely 
track the steep rise in magnetic field, with half the electron heating taking place on a scale 
of 17.3 km, corresponding to 6.4 electron inertial lengths and a small fraction (0.15) of an 
ion inertial length. Although much of the electron dynamics is linked to the DC electric and 
magnetic fields within the ramp (Feldma n et al.|198 3 ; Goodrich and Scudder 1984; Scudder 
1 1995 ; Lefebvre et al. 2007]) and is therefore reversible (the distribution function in this lim- 
iting case might be dependent upon energy and adiabatic invariant in de Hoffmann-Teller 
reference frame), the fact that both T e \\ and T e j_ rise together suggests an inflation of the 
particle phase space distribution that is not reversible, due primarily to the filling in and/or 
entrapment of electrons in regions of phase space that would otherwise be inaccessible. 

This infilling can be seen in the cuts of the distributions shown in Figure [22] Within the 
steep ramp, the inflated distribution is evident, with the flat-topped infilled region already at 
its downstream level. This supports the notion that the temperature profiles shown in Fig- 
ure [2T] really do represent irreversible heating. Interestingly, Figure [22] shows that features 
previously reported with the ramp, e.g., the beam vestige of the solar wind peak (Feldman 
et al. 1983), are present only in the more gradual initial rise that precedes the steep ramp. 
That beam has been totally eroded by the time this electron scale ramp is encountered. 

Thus the electron heating occurs over scales that are significantly smaller than the con- 
verted proton gyro-scale V n / '£l c % invoked in Bale et al. ( 2003 I and also smaller than the ion 
inertial length that might be anticipated due to micro-instabilities within the shock current 
layer ( |Papadopoulos| 1985b] |Matsukiyo and Schoier|2 006 1. 

Recent statistical studies (Hobara et al. 2010) argued that previous fits to a proxy of the 
plasma density profile ( |Bale et al.|2003| > mixed contributions from the more extended foot 
region governed by reflected gyrating ions. Restricting the measurements to just the steep 
ramp, they report widths in the range 3-55 c/ui pe with a decreasing trend as the Mach num- 
ber increases. They interpreted their work in terms of shock steepening limited by the dis- 
persion of electron whistler waves, with dispersion relation u) = Q C e cos#b„ (fc 2 c 2 /ajp e ). 
The limiting case of a wave capable of phase standing in the incident flow has a wavelength 
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Fig. 22 Cuts of the electron distribution functions in the solar wind, initial ramp, steep ramp, and downstream 
along (solid) and perpendicular (dashed) to the magnetic field. The locations of the cuts are indicated along 
the axes in Figs. |20| and |21| Note the solar wind halo drift evident in the anti-aligned direction and the absence 
of features within the steep ramp. (Adapted from Schwartz et al. 201 1 J. 



Table 1 Shock Parameters 2005 Jan 9 @ 22:15 



Parameter 


Value 




^shock 


+10.8 


kms 1 


Unshocked magnetic field B«' 


(3.07, 1.35, 8.14) 


nT 


Unshocked electron density 


4.0 


cm -3 


Location (Earth radii) 


(12.3, 13.3,-6.7) 


R e 


n shock normal (timing) 


(0.855, 0.418, -0.307) 




n (model)|Schwartz||l998} 
V n = V ■ n (shock rest frame) 


(0.904, 0.383,-0.189) 




373 


kms -1 


Alfven Mach no. Ma 


3.8 




Magnetosonic Mach no. M ms 


3.0 




9 Bnu = ZB, n 


83 


o 


Plasma ion ft 


0.4 




Plasma electron /3 e 


0.34 




Electron inertial length c/io pe 


2.7 


km 


Ion inertial length c/uj pi 


117 


km 


V n /Q clu » 


443 


km 




139 


km 


Whistler wavelength A 


24.8 


km 


Electron Larmor radius r^ e u 


1.01 


km 



t All vectors are in the GSE frame of reference. Subscripts "s" ("«") denote quantities in the (un)shocked 
region. 

tt Q ci = eB/m,p is the proton gyrofrequency 



that can be written 

A cos6) s „ [rrH 

= ^ ^7 \ I 

C/Wpe Ma V me 

The results from Table [T] yield a value of 9.2 for this ratio, comparable to the 6.4 electron 
inertial lengths given above. The fact that supercritical shocks steepen to this whistler limit 
suggests that dissipation processes are insufficient to broaden the transition further. 
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It should come as no surprise that the steepening of a fast mode (right-handed) wave 
results in a right-handed whistler signature. Indeed, the non-coplanar component of the 
magnetic field (Thomsen et al. 1987 ), responsible for the difference in the shock electro- 
static potential when viewed in different shock rest frames I Goodrich and Scudder 1984}, is 



right-handed. There is new evidence fSundkvist et al.|2012 1 that the wave Poynting flux is 



directed away from the ramp region upstream as expected for dispersion-limited steepening. 

The present study measures directly the actual temperature profile of the electrons. The 
result confirms that nonlinear steepening proceeds down to scales limited by whistler dis- 
persion. We have argued that this represents irreversible heating, implying that dissipation 
is operative on this, or probably smaller, scales. 

We have attempted a similar analysis on other shock crossings observed by Cluster, 
with consistent findings. Suitable events are rare, since they require the combination of a 
slowly moving shock and favorable magnetic field orientations. Future space missions need 
to be proposed to target electron physics and hence should provide numerous examples for 
statistical studies. 

What process(es) are actually responsible for (sub-)whistler-scale dissipation? The over- 
all inflation in phase space is linked to the action of the cross-shock electrostatic potential in 
concert with the magnetic mirror forces. Some or all of the potential may be concentrated 
in intense spikes (Bale and Mozer 2007) that may break the adiabaticity of electron phase 
space trajectories despite a ramp thickness which, in our example, is 20 times the local elec- 
tron gyroradius. It is worth noting that the localized spikes of the electric field are present 
inside the ramp region. Figure [23] represents 10 second interval of electric field measure- 
ments inside the ramp by Cluster 2, from 22:15:30 to 22:15:40 that is relatively short time 
with respect to time of the shock crossing but corresponds to ramp region (courtesy of F. 
Mozer). One can clearly see quite intense bursts of the electric field having amplitudes as 
large as 20-30mVm~ 1 . These bursts are very similar to those reported in Section [2] Such 
electric field bursts can be one of the possible sources of electron heating and scattering. 

Another candidate processes (e.g. Balikhin and Gedalin 1994 ) responsible for in-filling 
regions of phase space, in some of which electrons are trapped, include wave scattering 
(|Scudder et al.|1986c[|Veltri and Zimbardo|1993| ) and demagnetization (Balikhin a nd Gedalin| 



1994); these will require further analysis and simulations. 

Our discovery of short scale electron heating has an important consequence for electron 
and ion acceleration. Gradient drift and surfing mechanisms are sensitive to the scale of the 
field transitions (Zank et al. 1996), becoming very efficient at scales comparable to those 
reported here. 



6 Can anomalous resistivity account for energy dissipation and electron heating. 

The major results reported in this Section were first published in Balikhi n" et al.| ([2005 1; 
[Walker et al.|p0"0"8) . 



6.1 Ion Sound Wave Packets at the Quasiperpendicular Shock Front. 

The data used in both articles presenting observations of short scale waves were collected 
by the EFW instrument on board Cluster satellites. EFW uses two pairs of spherical probes 
in the satellite spin plane situated on the ends of wire booms whose length is 44 m as shown 
in the left hand panel of Figure[24] Thus, the distance between probes adjacent/opposite to 
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Fig. 23 Electric field measured onboard Cluster 2 satellite from 22:15:31 to 22:15:41. Electric field bursty 
spikes having amplitudes of 20-40mVm — 1 having duration of the order of 0.1 sec are clearly seen during 
ramp crossing. Figure is provided by F. Mozer. 

one another is ~ 62/88m respectfully. Normally, the EFW instruments return the electric 
field calculated as the difference in probe potentials between probes 1 and 2 (E12) and 3 and 
4 (£34) with a sampling rate of either 25Hz (normal science mode) or 450Hz (burst science 
mode). The individual probe potentials are also available with a time resolution of 5Hz. In 
addition to these standard modes, there is a triggered internal burst mode. Using this mode, 
data for a short time period may be captured with a much higher sampling rate. The EFW 
data that has been analysed in this study consists of internal burst mode data comprising the 
four individual probe potentials sampled at 9kHz for periods of around 10 seconds. Since 
the internal burst data is captured and stored depending upon some criteria, it may be that 
although the shock region was targeted for data collection, the waveforms returned may not 
have been captured in the shock front itself. To this end, a search was made to find possible 
candidate events by cross referencing the list of Cluster shock crossings for 2002 with the 
list of periods for which internal burst data are available. This resulted in a list of 10 possible 
events. Of these events, a comparison between the FGM magnetic field measurements and 
the time periods for which EFW internal burst data were available showed that there were 
only two shocks for which the period of internal burst data lay solely in the foot region of 
the shock. Of these, one shock possessed a magnetic profile that was highly turbulent and 
difficult to interpret and was also eliminated from further analysis. This left just one clean 
shock on which to perform the analysis. 

The internal burst data sets are the only ones generated by EFW that contain the indi- 
vidual probe potentials at a high enough sampling rate to investigate waves and turbulence 
at frequencies around the lower-hybrid frequency in the vicinity of the terrestrial bow shock 
(10-30Hz). By using the individual probe potentials it is possible to compute two parallel 
electric field components one on either side of the satellite. For example, the probe pairs 1,3 
and 4,2 maybe used to compute electric field components S13 and £42 whose directions are 
parallel and are spatially separated by a distance of ~62.2m. This technique has previously 
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Fig. 24 The left panel shows the configuration of the EFW electric field probes and illustrates the electric 
fields calculated from them. The right hand panel shows the magnitude of the magnetic field (bottom) and 
ion bulk flow (top) measured during the shock crossing that occurred on February 26 th , 2002 at 2134 UT. 
(Adapted from Balikhin et al.|2005) . 



been used by Balikhin et al. (2005} and |Tjulin et al.| ( |2003) > to study small scale electric 
field structures and waves and is similar to the short baseline interferometry techniques em- 
ployed in the analysis of data from sounding rockets (Pinc onet al.|1997) . Since the probe 
potentials can be used to calculate two parallel electric field vectors it is possible to use the 
phase differencing technique to determine the wave vector k. This method may also be used 
to examine the polarisation characteristics of the wave in question. In this case, the phase 
differencing algorithm is applied to a pair of perpendicular components of the electric field 
(as opposed to the parallel field components mentioned above). The resulting histogram of 
the phase difference as a function of frequency yields a vertical line of constant phase dif- 
ference with respect to frequency at a phase difference of zero for a linearly polarised wave 
and ±7r/2 for a circularly polarised wave. Thus, this technique may be used to help dis- 
tinguish between a linearly polarised lower hybrid wave and a circularly polarised whistler 
mode wave, both of which have been observed at these frequencies. This method is used in 
preference to an examination of the coherency (see for example Krasnosels kikh et al.|1991| > 
due to the short duration of the wave packets. 

The magnetometer data, used to put the electric field measurements into context within 
the shock front and compute the lower hybrid resonance frequency, come from the FGM 
instruments (Balog h et al.|1997] > and made publically available through the Cluster Active 
Archive. These measurements typically have a sampling rate of 22Hz. 

All the data presented in this Section were recorded during one shock cros sing on Febru- 



ary 26 i/l ,2002 at around 2134 UT during the time intervals marked on Figure 125] by vertical 



lines. Red lines mark the periods of registration of ion sound waves, green lines the periods 
of registration of lower hybrid and whistler waves. 

During this period the Cluster satellites were situated in the foot region of a quasiper- 
pendicular shock {Bbu ~ 55°, Ma ~ 4.3). The EFW instrument onboard Cluster 3 was 
triggered to operate in internal burst mode for a few seconds. 

Two parallel electric field vectors of these electric field measurements lie in the same 
direction and have a perpendicular separation of ~ 62.2m in the direction Pi to P3. The 
availability of two closely spaced, simultaneous measurements enables the use of phase dif- 
ferencing techniques (Balikhin et a l.|1997a| l for the identification of propagation modes for 
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Fig. 25 Waveforms of electric field measurements during shock front crossing on February 26 th , 2002 at 
2134 UT. The vertical lines mark periods where the waves were registered. As it will be shown later red 
vertical lines mark periods where the waves were identified as ion-sound, two first green columns as lower- 
hybrid electrostatic waves, third green column as whistler waves in lower hybrid frequency range. 

waves with coherence lengths down to a few Debye lengths based upon single satellite mea- 
surements. Since there is no component measured normal to the spin plane, the separation 
between temporal and spatial variations is possible only in the spacecraft spin plane. As a 
consequence, phase differencing methods are limited to the determination of the projection 
of the k-vector in the spin plane. In most cases, however, this can provide enough infor- 
mation to identify the plasma wave mode. This approach was implemented in these studies. 
Plasma measurements were obtained from the CIS HIA (ions) and PEACE (electron) instru- 
ments. Magnetic field data were obtained from FGM. It should be noted that the spin vector 
of the Cluster satellites is almost coincident (to within 5°) with the z GSE axis. 

The ion bulk velocity (top panel) and the magnitude of the magnetic field (lower panel) 
as measured by Cluster 3 spacecraft are plotted in the right hand panel of Figure[24] Initially, 
the spacecraft was in the solar wind. The foot region was encountered just before 2134UT 
and the shock ramp was crossed around 2134:1 2. 5UT. The plasma bulk velocity began to 
decrease around 2133:50. Shortly before 2134 low frequency oscillations were observed in 
the magnetic field, a feature commonly observed in the foot region of supercritical shocks. 
The beginning of the foot region is characterised by a large amplitude, nonlinear structure 
similar to those previously reported by Walker et al. ( 1 999a ). A comparison of magnetic field 
and plasma data show that this structure is not a partial penetration of the ramp. The present 
study is limited to the short interval at the beginning of the internal burst mode indicated by 
the vertical line and coincides with the foot region. 

The electric field component £31 as measured during the initial part of the internal burst 
mode interval is shown in the upper left panel of Figure [26] and its Morlet wavelet spectra 
is shown in the lower left panel. The electric field fluctuations show a pair of well defined 
wave packets centered around 2134:01.6 and 2134:02.05UT. Their frequency ranges are 
100-800Hz and 250- 2000Hz respectively. We present here the results of the identification of 
these wave packets to illustrate the use of the technique and its results. 

The / — &23 spectrum, as shown in the center and right hand panels of Figure [26] is a 
histogram representation of the distribution of wave energy in frequency-fc space for the first 
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Fig. 26 Left: The waveform (top panel) and wavelet spectrogram (bottom panel)of the electric field computed 
from the difference in potential between probes 3 and 1 . Centre and right: Examples of the f — k spectrograms 
for the first (left) and second (right) wave packets. (Adapted from Balikhi n et al.|2005) 



Fig. 27 A comparison of the wave vector directions for the two wave packets. The dotted and dashed lines 
represent the X component of wave vector for events 1 and 2 respectively. The corresponding Y components 
are shown by the dash-dotted and solid lines. (Adapted from Balikhin et al. 2005) 



wave packet (Balik hin et al.|199 7al. The f — k spectrum shows a well developed ridge like 
maxima, the shape of which indicates the wave dispersion relation projected along the A?23 
direction. This result may be combined with a similar dispersion along the k±z to yield the 
wave vector projection in the satellite spin plane. 

Since the angle between the spacecraft spin plane and the GSE XY plane is small, we 
will consider that the projection into the spin plane is the same as that into the GSE XY 
plane. The projection of the dispersion relation into the GSE XY plane is shown as the 
solid line in Figure [27] for the first (centre anel in Figure |26|l and second (right on previous 
Figure) wave packets. The observed frequency range of the first wave packet (100 — 800Hz) 
corresponds to approximately 0.25 — 1.9J7 ce , and the magnitude of wave vector projection 
is in the range 0.015 < fci < 0.075m _1 . For this interval the electron temperature is T e ~ 
17eV and plasma density = 9.7cm -3 . This leads to an estimate for the Debye length of 
\ d w 10m. Thus the observed values of for the projection of k correspond to ^ 8 — 40A d . 

The satellite frame dispersion relation in the satellite spin plane is shown by the solid 



line in Figure 27 It's phase velocity is in the range 40-70kms~\ The Doppler shift can 
be estimated as the scalar product of the solar wind velocity and spin plane wave vector 
component. This estimation of the Doppler shift term is shown as a dashed line. It has the 
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same sign as the phase velocity and is always greater than the observed wave dispersion 
indicating that in the plasma frame the waves propagate in the direction opposite to that of 
the solar wind, but are convected Earthwards by the plasma flow. This convection reverses 
the direction of propagation in the satellite frame. The average angle between the spin plane 
projections of wave vector and the magnetic field is about 20°. 

The second wave packet analysed was observed « 0.3 seconds after the first. The elec- 
tric field waveforms (not shown) again indicate a good correlation between the correspond- 
ing electric field components measured by different probe pairs. The / — /C23 spectrum 
calculated for this wave packet is shown in the right panel of Figure [26] The ridge like 
maxima in these spectra correspond to the projections of the wave dispersion relation in 
the direction ^23. The resulting dispersion relation is shown as the solid line in Figure [27] 
. Its frequency range is 250 — 2000Hz (w 0.6 — 4.9J? ce ), and the magnitude of wave vec- 
tor projections is in the range « 0.018 - 0.075m" 1 . For this wave packet, the satellite 
frame phase velocity is in the range ISO-ieOkms" 1 . The range of wave vectors and angle 
of propagation with respect to the magnetic field for the second wave packet coincide with 
those determined for the first. Even more surprising is the fact that the angle between the 
two wave vector projections is less than 5° . The dashed line in Figure 27 shows the esti- 
mation of the Doppler shift. It can be seen that the Doppler shift term for the second wave 
packet is less than that of the observed frequency and so the second wave packet propa- 
gates in the same direction in both the satellite and plasma frames. Therefore the first and 
second wave packets propagate in opposite directions in the plasma frame. While for the 
second wave packet the satellite frame phase velocity is the sum of of its plasma frame ve- 
locity and the solar wind convection speed for the first wave packet it is their difference. 
That explains why in the satellite frame the second wave packet propagates faster than the 
first one. The use of multipoint measurements enables one to separate temporal and spatial 
variations. In the current study it is possible to distinguish which of these two wave modes 
was observed. Thus we have a method that is independent of using the observed frequency 
criterion formulated by Gurnett l |1985| >. For this interval \B\ ~14.8nT and hence the lo- 
cal electron cyclotron frequency f e = Q C e/2pi ~415Hz. As can be seen from the / — k 
spectra shown in Figure [26] that the maximum wave energy of the first wave packet occurs 
at a frequency lower than f e . According to the classification used by Gurnett ( 1985) this 
should be a whistler wave packet whose dispersion relation may be written as (neglecting 
thermal corrections) J 2 = Q 2 e cos 2 6 Bk k 2 c 2 /(k 2 c 2 +cjp e ), where 9 Bk , Q e , the 
angle between the wave vector and the magnetic field, the electron cyclotron and electron 
plasma frequencies respectively. The wave vectors for the first wave packet lie in the range 
(kc/ui) 2 ~ 30—150 and therefore correspond to the electrostatic limit of the mode for which 
the plasma frame frequency should be ~ co e cos 8 ' B k- If we estimate the angle 9sk using 
the angle between the projections of wave vector and the magnetic field in the spin plane, 
the plasma frame frequency can be estimated as f e cosdsk ~280Hz. For the wavevectors 
found, 0.015 < k < 0.0075m~ 1 the electrostatic whistler phase velocity varies in the range 
24 < v p h < 112kms _1 in the plasma rest frame. In the spacecraft frame the slowest waves 
would reverse their direction of propagation, so that waves propagating in both directions 
would be observed. However, it has been shown earlier that all waves are propagating in the 
same direction. Moreover, for the strongly dispersive electrostatic whistler the phase veloc- 
ity should vary by a factor of two or more over the observed range of wavevectors, while 
the actual variation is within 20% only. These arguments exclude the possibility that the 
observed mode is the whistler in the electrostatic regime. 

The other possibility is the ion-sound mode. Since we are limited to spin plane mea- 
surements of wave vectors only order of magnitude estimations of the wave parameters 
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can be made. For such crude calculations we will disregard the factor Ogk ~ 18° in dis- 
persion of ion-sound waves and use the simplified form lu = kvi s / Jl + k 2 \ 2 D where 

vis = \/kbT e /mi is the ion-sound velocity and k b is Boltzmann's constant. During the time 
interval in which both waves packets were observed v is « 40kms _1 . Thus, in the plasma 
rest frame the wave phase velocity should be in the range 0.80?;^ < v p h.pf < 0.99«i s . This 
velocity dispersion is very close to the observations. If observed waves are indeed ion-sound 
waves their plasma frame frequency should be in the range « 75—100 Hz, much lower than 
the observed frequency. This disagreement can be attributed to the Doppler shifts estimated 
as ^ Vsw ~ 600 — 3000Hz. In reality the Doppler shift is smaller due to the angle between 
the wave vector and the solar wind velocity. 

The above arguments indicate that the first wave packet consists of ion-sound waves. As 
previously mentioned, the wave vectors for the second packet have exactly the same range 
as the first. Therefore, all arguments used above to deduce the wave mode of the first wave 
packet are valid for the second. The main difference between these two wave packets is in 
the sign of the Doppler shift. For the first wave packet, the observed frequency is the differ- 
ence between the Doppler shift and the plasma frame frequency whilst for the second it is 
their sum. It can be seen that they almost coincide for the whole range of observed waves. 
The angle between the averaged propagation directions of these wave packets is < 5°. This 
coincidence in the parameters for these two wave packets, observed at clearly distinct pe- 
riods of time can only be explained by their simultaneous generation at the same location. 
The generation of ion-sound waves at the shock front are usually attributed either to electric 
currents or the strong electron temperature gradients in the ramp. Both waves packets were 
observed upstream of the ramp and carried by the solar wind flow towards it. Since there 
appear to be no strong gradients in the electron temperature in the foot these waves are prob- 
ably generated by electric currents. The very short duration of these waves indicates that the 
current layer might be localized in space and time. Such small scale current layers have been 
predicted by a nonstationary model of the shock front ( |Krasnoselskikh|1985| |Galeev et al.| 
[T988l[T989l|Balikhin et al.|1997b||Walker et al.|1999a| l. In this model quasiperiodic steep- 
ening of and overturning of the shock ramp takes place leading to the ejection of a nonlinear 
whistler wave into the upstream region. The amplitude of these nonlinear structures can be 



comparable to the |B| changes in the ramp itself (Walker et al. 1999a| > and will be associated 



with localised currents responsible for the ion-sound waves. 



6.2 Observations of lower-hybrid waves 

The data set used in this study was collected by the EFW instrument in the same burst mode 
regime as in previous case onboard the Cluster satellites using onboard timing provided by 
the DWP instrument l |Woolliscroft et al.|1997) . 

Figure |28fl) shows an overview of the magnetic profile of the shock encountered on 



February 2fr^ at 21:34UT. From Figure 28 1) it can be seen that Cluster 3 first encountered 
the foot region of the shock just before 21:34UT, finally crossing the ramp and entering the 
downstream region at approximately 2134: 15UT. Here the EFW internal burst data selection 
was triggered at 21:34:01.922UT and lasted for a period of 10.47 seconds as indicated by 
the shaded region in the figure. 

The analysis presented here was performed on data recorded on February 26 th , 2002 
just after 21:34UT on spacecraft 3. This quasi-perpendicular shock crossing took place on 
an inbound pass at a position (12.0, -1.60, 8.07)Rg. As can be seen from Figure |28f 1) 
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Fig. 28 Panel(l): The magnetic profile of the bow shock crossing observed by Cluster 3 on February 26 th , 
2002 just after 2134UT. The period for which EFW internal burst mode data is available is indicated by the 
shaded region. Panel(2): The wavelet dynamic spectrogram of electric fields E12 (top) and E34 showing a the 
occurrence of oscillations just above the lower-hybrid resonance frequency (black line) for event 1 . (Adapted 
from |Wa!ker et al.|2008> . 



the whole period of internal burst data was collected in the foot region of the shock. The 
initial increase in the magnetic field profile at around 21:34UT as has been shown above and 
published in Balikhin et al. (2005) to be part of the foot region rather than a partial ramp 
crossing. 

During the 10.5 second period for which EFW internal burst data is available there 
were several short periods when the electric field measurements indicated that there were 
oscillations occurring at or just above the local lower-hybrid resonance frequency. In the 
following subsections the properties of the waves observed are discussed. 

The first event occurred just after 2134:05UT. Figure [28f 2) shows a dynamic spectro- 
gram of the electric fields E12 (top) and E34 measured between probes Pi and P2 and probes 
P3 and P4 respectively calculated using a Morlet wavelet transform. The black line repre- 
sents the lower-hybrid resonance frequency. It is clearly seen that at around 2134:05.2 and 
there is a packet of waves at a frequency between 10-20Hz, whose lower edge is just above 
the lower-hybrid resonance frequency. The duration of this wave packet is around 3ms cor- 
responding to a few wave periods. Having identified a possible occurrence of lower hybrid 
waves, the phase differencing technique was applied to parallel electric field vectors in an 
attempt to compute the dispersion relation of the waves and hence provide an unambiguous 
identification of the wave mode. However in this case using the spin plane electric field com- 
ponents £13 and Ei2 in the frequency range of interest (10-20Hz) no measurable dispersion 
of the waves was observed on scales of the separation distance of the probe pairs (62.2m), 
see Figure |29|l). 

This implies that the wave travels over the spacecraft at rather high speed so that there 
is virtually no difference in the phase of the wave measured at the two points on either 
side of the satellite. This was also evident in the waveform of the electric field signals. A 
comparison of the waveforms (Figure [29{2)) shows that the two measurements which are 
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Fig. 29 Panel(l): The u> — k joint spectrum calculated from the phase differences measured between the two 
parallel electric field components E13 and E42 for event 1. Panel (2): The electric field waveforms E13 and 
E42 (top panel) and E41 and E23 (lower panel) for the first event. (Adapted from Walker et al. 2008). 

observed to vary in phase which indicates that whatever passed over the satellite has a scale 
much larger than the individual probe separation distances. 

So, one should conclude that the phase differencing method appears to be unable to 
show any dispersion in the parallel electric field vectors that means that this method cannot 
be used to reliably identify the wave packet as being lower hybrid. As a result, one needs 
to investigate some other wave properties of the wave packet to determine if they are com- 
patible with the lower hybrid mode. As it was mentioned above that by applying the phase 
differencing method to perpendicular components of the electric field it should be possi- 
ble to determine whether the wave packet is linearly or circularly polarised. To this end 
one should calculate the phase difference between two pairs of probes. The result of this 
calculation gives the estimate that in the frequency range 10-20Hz the value of the phase 
difference is around zero. This result indicates that the wave possesses linear polarisation. 
This leads to the conclusion that the wave mode that is observed in this case corresponds to 
a lower-hybrid mode. 

Similar analysis of the second event occurred around 2134:04.5 on February 26 th , 2002 
shown the very same result, namely, it shows the phase difference of zero which again 
indicates that the wave is propagating with a large phase speed over the satellite and the 
wave packet possesses a linear polarisation and it exhibits properties that are consistent with 
propagation in in the lower hybrid mode. 

The third event highlighted by the authors of the paper | Walker et al.| p008[ l occurred 
between 2134:07.3 and 2134:07.45 UT on February 26 th , 2002. The wavelet dynamic spec- 
trogram analysis showed a wave packet in the frequency range 10-15Hz that lies just above 
the lower hybrid resonance frequency. This wave packet was observed to drift in frequency 
as time increases. This change in frequency mirrors the change in the lower hybrid reso- 
nance frequency as calculated from the magnetic field. Analysis of parallel electric field 
components using the same phase differencing method again indicated zero phase differ- 
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ence between the components. However, a comparison of perpendicular electric field com- 
ponents indicated that a phase difference between two signals is equal to tt/2. This leads to 
conclusion that the wave packet possesses circular polarisation and is thus not propagating 
in the lower hybrid mode. The circular polarisation indicates that this particular wave packet 
is propagating in the whistler mode. 



6.3 Estimates of efficient collision frequency using direct measurements of ion-sound and 
lower-hybrid waves 

A definition of the problem of conductivity relies on exchange of momentum between elec- 
trons and waves assuming the current is mainly carried by electrons. The conventional for- 
mula for plasma conductivity reads 

ne 2 
m e v 

where n is the plasma number density and v is the collision frequency of electrons with 
scattering centers, usually ions or neutrals with respect to momentum loss. When electrons 
excite some oscillations or waves as a result of instability development they also loose the 
momentum and this loss is referred as the anomalous momentim loss. In order to find effec- 
tive collision frequency v e ff one has to use the momentum conservation law in the system 
consisting of electrons and waves. In the case of instability this momentum exchange can be 
written as follows 



v e ffrn e n lt ed = ^ J d k 1kW k ( — ) 
where Tt e( i is the relative velocity of electrons carrying current, j k is the instability incre- 

\E\ 2 

ment, W k is the wave energy density that is defined as W k = °' 2 1 , | E | is the turbulent 
electric field amplitude. We have in the left hand side the rate of the electron momentum 
loss per unit time, and in the right hand side we have the momentum gain by waves due to 
instability. It follows then that 



v e ff = 775-^3 n=5 — r / d k-y ke W k { — ) 

(2-K) 6 m e n \ Tt ed | J cu k 

If one would like to evaluate the efficient collision frequency from direct measurements 
it is necessary to have an estimate of the wave energy. Using this estimate one can evaluate 
the shock thickness that shock might have if it would be determined by the anomalous 
collisions using the characteristic length of the momentum loss: 

Vsw 



"iff 



The comparison of the thickness obtained from this estimate with the real shock thick- 
ness can be used to evaluate the relative role of the efficient collisions. 

The standard estimate of the efficient collision frequency for ion sound mode ( |Galeev| 
|and Sagdeev|1984| reads 



W 



' noksTe 

Taking the estimate of the averaged electric field intensity (E) 2 » 1CP 5 calculated mak- 
ing use the data of measurements (it varies from %1 x 10 to 2 x 10~ 5 V/m) and density 
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and temperature from observations n = 9.7cm -3 , T e = 17eV, oj p = 1.7 x 10 5 sec -1 , one 
can find W = 4.5 x 10~ 17 ,W/{nk B T e ) = 1.8 x 1CT 6 , 

ul%f = 0.3s _1 ,L a n = ^ SW ■ ~ 1200fcm 
V eff 

, where v 1 ^ f is anomalous collision frequency due to ion-sound wave activity, L an is the 
characteristic scale of anomalous energy exchange between electrons and ions. It is suffi- 
ciently larger than the electron inertial scale c/uj pe = 1.76km, and comparable with the 
thickness of the foot region. 

Another group of waves, namely lower hybrid has maximum amplitudes of the order of 
10mVm -1 and average electric field energy density of the same order of magnitude as ion 
sound waves (from 1 x 1CF 5 to 4 xlO~ 5 V/m). In order to evaluate the efficient collision 
frequency for these waves one should take into account the properties of lower-hybrid drift 
waves. To this end we shall rely on the study of lower hybrid drift instability published by 
|Davidson et al.| ( [T977| l. The maximum growth rate for these waves can be estimated as 

1LH ~ OL^LH 

where Qlh is lower hybrid frequency, coefficient a < 1, typically a ~ 0.1, and can reach 
values up to 0.3. Taking maximum of the linear growth rate we can evaluate the upper limit 
of the effective collision frequency. The phase velocity of waves around the maximum of 
increment is of the order of ion thermal velocity of ions, and the drift velocity of electrons 
that carry the current can be estimated evaluating current velocity from macroscopic gradient 
of the magnetic field. This estimation gives the value comparable with ion thermal velocity. 
Thus the estimate of the efficient collision frequency in this case can be written as follows: 

M anLH m w 

ejtj * m e nk B Ti 

In the region of observations where B = 14nT and lower hybrid frequency is approxi- 
mately equal to 56sec~ 1 , thus the efficient collision frequency for these waves is found to 
be of the of the order of 

v eff w q^l/j— ~ 0.1 x 2 x 10 3 x 56 x 10~ 6 ~ 10 _2 s _1 

that is sufficiently smaller than the efficient collision frequency for ion sound waves. 
The characteristic dissipation scale 

Lan ~ TTT-n km w 3.5 x 10 4 fcm 
10 z 

and is sufficiently larger than the major characteristic scales of the shock front. 

This leads to the conclusion that the anomalous resistivity observed can not account 
for the important dissipation rate. The characteristic scales of the dissipation are too large 
compared to the shock transition features observed. 



7 Nonstationarity and reformation of high Mach number quasiperpendicular shocks: 
Cluster observations 

In this Section, using Cluster observations, we provide convincing evidence that high-Mach- 
number quasiperpendicular shocks are indeed nonstationary, and moreover, quasi-periodic 
shock front reformation takes place. Most of the material of this Section was first published 
in |Lobzin et^ ( |2007l l. 
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7.1 An Example of a Typical Crossing of Nonstationary Quasiperpendicular Shock Wave 

A number of magnetic field profiles of the quasiperpendicular terrestrial bow shock ob- 
served by Cluster triaxial flux gate magnetometers (FGM) (Balog h et al. 11997] ) in the period 
January-May 2001 were studied. It was found that nonstationarity seems to be typical for 
shocks with relatively high Mach numbers. Both from numerical simulations and exper- 
iments it follows that the details of this nonstationary behaviour of the shock front may 
depend strongly not only on the fast magnetosonic Mach number, Mf, but also on the up- 
stream /3 6) j and the angle between the upstream magnetic field and the shock normal, #b„. 
For a detailed case study, a shock was chosen that could be considered as a typical quasiper- 
pendicular, supercritical, high-Mach-number shock wave, namely the shock crossing that 
occurred on 24 January 2001 at 07:05:00-07:09:00. Indeed, from the available experimental 
data and with the use of the multi-spacecraft timing algorithm described by Schwartz (1998 1 
the following estimates were obtained: /3 e = 1.7, ft = 2.0, 0B n = 81°, Ma = 10, and 
M f =5. 

Figure [30] shows the magnetic field profiles measured by the Cluster FGM instruments 
on January 24 200l"\ 2001. The panels on the left show the full resolution data, sampled 
at 67 Hz (black line) and the result of averaging this data using a 4 second sliding window 
(red line). The panels on the right show the result of low pass filtering the data at 2Hz. 
This process enhances any large peaks in the magnetic field measurements. All the profiles 
can be considered as quite typical for high-Mach-number quasiperpendicular shock waves. 
From the averaged data shown by the red lines we observe that the shock front consists of 
a foot, a ramp, and at least one overshoot-undershoot cycle, i.e. large amplitude peak of 
the magnetic field at the end of the ramp region and following after it minimum. The small- 
scale oscillations of large amplitude are superimposed on this large-scale structure. To check 
whether these fluctuations are consistent with plane wave activity, the degree of polarization 
for the magnetic field waveforms obtained from STAFF experiment ( Cornilleau- Wehrlin 
|et al.| 1997] >. By definition, the degree of polarization approaches a unity if and only if most 
of the energy is associated with a plane wave (Samson and Olson 1980). It was found that 
between the forward edge of the shock and the magnetic overshoot the oscillations in the 
frequency range 3-8 Hz have a high degree of polarization greater than 0.7 and that this po- 
larization is elliptical. This wave activity can be considered as a whistler wave train nested 
within the shock l |Galeev et al.|1988[[T9g9)|Krasnoselskikh et al.|2002| ). Obviously, the pres- 
ence of whistler oscillations, due to their high amplitude, has a considerable impact on the 
large-scale shock structure. Indeed, averaging of magnetic field data reveals two regions re- 
sembling overshoots for SC4 whilst only one maximum is observed for SCI. The profiles 
for the other spacecraft appear to be more complicated. It follows from these considerations 
that the concepts of both overshoot and ramp, which must precede it, become ambiguous for 
such nonstationary shocks. Instead, we can speak about short scale large-amplitude struc- 
tures embedded into the shock transition, with the forward edge of one of these structures 
playing a role of the ramp. Figure [30] also shows that the magnetic field profiles measured 
onboard the different spacecraft differ considerably from each other. Obviously, the number 
of large-amplitude peaks, their amplitudes, as well the positions within the shock front, are 
different. The waves observed by different spacecraft in the foot region are also different. In 
particular, from Figure[30](left) it is easily seen that the time interval between the beginning 
of the wave activity at the forward edge of the shock and the ramp crossing may differ by 
10-20 s. This difference is substantial compared to the duration of the crossing of typical 
elements of the shock structure. The distinctions found between observations from the dif- 
ferent spacecraft are related to temporal rather than spatial variations in the structure of the 
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Fig. 30 The magnetic field profiles obtained by FGM experiments aboard four Cluster spacecraft during the 
Earth's bow shock crossing on 24 January 2001. (left) High-resolution magnetic field data (black line) and 
the data obtained by sliding averaging over 4 s time intervals (red line), (right) Vicinity of overshoots, with 
large peaks in the magnetic field magnitude. Oscillations with frequencies higher than 2 Hz were removed. 
To emphasize the similarity and differences of the profiles, the data for the first 3 spacecraft are shifted with 
respect to that for the 4th one. (Adapted from Lobzin et al. 20071. 



shock front because the spacecraft separation is comparable with shock front thickness. In- 
deed, the distances between spacecraft lie within the range 380980 km. The foot thickness 
estimated with the use of the theoretical formula derived by Sch wartz et al.| ( |T983[ ) is equal 
to 550 km, in reasonable agreement with the observations, while the total shock front thick- 
ness is considerably larger. On the other hand, the maximum time lag between the crossings 
is about 3Tbi, where T^i is the ion gyroperiod T^i = 15.5 s. This time lag is larger than 
the period of the shock reformation. Relying on theoretical considerations and results of 
numerical simulations, Krasnosels kikh et al.| ( |2002) l argue that this type of nonstationarity 
is closely related to nonlinear whistler wave trains embedded into the shock front and that 
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Fig. 31 Hodograms of magnetic field around their maxima for different satellites in the vicinity of supposed 
ramps. 



this is a typical property of quasiperpendicular high-Mach-number shocks. Further evidence 
for the existence of whistler waves embedded within the shock front can be seen from the 
rotational features of the magnetic field observed in the vicinity of the peaks, as shown in 
Figure [31] that are typical of whistler mode waves. 

The large-amplitude structures seen in the magnetic field profiles within the overshoot 
region and its vicinity have a characteristic time of about 2 s. To examine both the similar- 
ities and differences of these profiles, oscillations with frequencies higher than 2 Hz were 
removed by low pass filtering the data. The filtered data were then used to calculate a set 
of optimal cross-correlation coefficients for profile fragments that last 35 s and include a 
portion of foot and the entire overshoot region. The highest correlation was found between 
SCI and SC2, while the lowest one was between SC3 and SC4, a result that is in accordance 
with visual observations of the shifted profiles shown in Figure [30] (right). An additional 
analysis of the relative position of the spacecraft tetrahedron and the shock reveals that the 
similarity of the shock profiles seems to depend mainly on the time interval between the 
shock crossings and/or the spacecraft separation measured along the shock normal rather 
than on the distance along the shock surface which is in accordance with the interpretation 
that the observed variations are temporal rather than spatial. 

Further evidence favoring the nonstationarity of this bow shock crossing comes from 
WHISPER measurements. In passive mode this experiment provides electric field spectra of 
natural emissions in the frequency range 2-80 kHz ( Dec reau et al.|1 997). The frequency-time 
spectrogram obtained by WHISPER experiment aboard SCI is shown in Figure[32] together 
with the magnetic field profile with the same time scale. The bow shock crossing can be 
identified by a substantial enhancement in the electric field fluctuations within the frequency 
range 25 kHz. For SCI, maximum intensity for these oscillations is observed at 07:06:48 
UT. One of the most obvious features of these spectra is the presence of intense waves in the 
vicinity of the plasma frequency, f pe = 27 kHz, together with downshifted oscillations. The 
most intense feature is a narrow-band Langmuir emission with a frequency in the vicinity 
of f pe . As compared with Langmuir waves, the power density of downshifted oscillations 



0.4 



-0.4 



of -<^__^>, CO 



0.4 

of 



O «: 



■ i 1 1 

-0.8 -0.4 o 0.4 0.8 



0.4 - 

- 

-0.4 - 

I i I i 

0.4 0.8 I 

B 2 



0.8 



0.4 



III 

0.4 0.8 



0.4 

m" o - 

-0.4 - 



~i — ' i — i — i 

-0.8 -0.4 



0.4 

crT o 

-0.4 



Dynamic quasiperpendicular shock 



51 



Mil 

In 




Fig. 32 (top) Electric field spectra and (bottom) magnetic field profile obtained during the Earth's bow shock 
crossing on 24 January 2001 aboard SCI. The frequency-time spectrogram is measured by the Whisper 
experiment. The vertical white bands correspond to the time intervals when no data were obtained in the 
natural wave mode. The wave intensity is colour coded with the reference level of 10 -7 V rms /Hz 1//2 , where 
rms is the root mean square to notify the averaged level of the fluctuating electric field variance. The magnetic 
field profile is obtained by FGM experiment. The time scales for the both panels are the same. (Adapted from 
|Lobzin et al.|2007) . 



is usually smaller, while the frequency band they occupy is considerably wider and can 
be as large as 15-20% of the central frequency. Both the plasma waves and downshifted 
oscillations are considered to be typical of the electron foreshock region. It is commonly 
believed that Langmuir waves are generated by a plasma-beam instability, while for the 
downshifted oscillations two different mechanisms have been proposed, namely, the plasma- 
beam interaction, see (see |Lacombe et aT] 1985; Fuse lier et al.||1985] l and the loss-cone 
instability of electron cyclotron modes ( [Lobzin et al.||2005| >. The mean frequency of the 
downshifted oscillations is not constant but varies within the range 0.2 - 1.0 f pe - In addition, 
there exists a tendency for a large shift to occur in the vicinity of the shock front, while 
near the edge of the electron foreshock the shifts are considerably smaller. However, this 
tendency exists only on large time scales of about 1.0-1.5 min. For smaller scales, 10- 
15 s, there are the large-amplitude variations of the mean frequency of the downshifted 
oscillations. The peculiarities of the spectra described above can be explained as follows. 
The downshifted oscillations are produced by energetic electrons, which are reflected by 
the bow shock and move almost along the magnetic field lines. Because the solar wind 
is quiet during the time interval considered (indeed, Figure [30] and Figure [5T| show that 
there are no significant variations of the magnetic field; the plasma bulk velocity is also 
approximately constant in the foreshock), the observed evolution of the wave spectra can 
only be attributed to variations of the suprathermal electron fluxes which are reflected from 
the bow shock and form the rabbit ears in the electron distributions upstream of the shock 
as was shown by |Lobzin et al.|(|2005|l. The reflecti on of electrons by a nearly perpendicular 
bow shock was studied by Leroy and Mangeney ( 1984) and |Wu| l |1984| >. They argued that 
the main characteristics of the distribution function of the reflected electrons depend first 
of all on the angle between the shock normal and upstream magnetic field, , and to a 
lesser extent on the ratio of the maximum magnetic field to its upstream value and on the 
electrostatic potential jump in the de Hoffmann-Teller frame. Resulting from shock front 
nonstationarity, slow variations of the effective normal of the reflecting part of the shock 
will lead to considerable variations of number density, energy of reflected electrons, and/or 
loss-cone angle, thereby producing the observed variations of the downshifted wave spectra. 
Both theoretical considerations and numerical modeling show that a characteristic time of 
the shock front oscillations or reformation is comparable with the ion gyroperiod (see |Leroy| 
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|et al.|1982[|Krasnoselskikh et al.|2002||Scholer et al.|2003[ >. The time scale of the spectra 
variations is also comparable with ion gyroperiod T^i, in accordance with our interpretation. 



7.2 Evidence for Shock Front Reformation 

As noted above, the magnetic field profiles for the shock under consideration have several 
nonstationary features. In this section, we consider large-amplitude structures, with a char- 
acteristic time of about 12 s and present the arguments in favor of front reformation for this 
particular bow shock crossing. Figure [30] (right) shows the magnetic field profiles obtained 
after low-pass filtering and shifting the data in time to clearly show the correspondence be- 
tween the elements observed aboard different spacecraft. For three spacecraft there are two 
large narrow peaks in the overshoot region and its vicinity, while for SC3 there is only one 
peak in the corresponding region (see Figure [30] (right), where these peaks are shown by 
arrows and numbered). The amplitudes of these peaks, both absolute and relative, differ for 
different spacecraft. In addition, the distance between two adjacent peaks also varies, be- 
ing the smallest for SC4 and the largest for SC2. Moreover, the single peak observed by 
SC3, which largest amplitude and relatively large width, may be formed due to the coales- 
cence of two separate peaks. The observed peaks in the overshoot region can be considered 
as a part of the nonstationary whistler wave packets since their rotational properties are 
clearly evident in Figure [3T| These properties were argued to be an intrinsic element of 
the quasiperpendicular supercritical shock front structure ( Krasnosels kikh et al.||2002| >. In 
order to investigate these features further, an analysis of their polarization was performed 
using the minimum variance technique. The results provide additional evidence in favor of 
shock front nonstationarity. Indeed, the corresponding elements have different hodograms, 
which can be rather complicated. However, some of the elements have approximately cir- 
cular polarization typical for large-amplitude whistlers as was stated in theoretical papers 
l |Galeev et alT][T988] [1989] |Krasnoselskikh et al.|[2002] l and is evidenced on Figure [3T] A 
comparison of the magnetic field profiles, shown in Figure 30 with the results of numerical 
simulations of high-Mach-number shock reformation (Krasnoselskikh et al. 2002) reveals a 
doubtless resemblance between them. Indeed, for large Mach numbers, quasiperiodic ref- 
ormation of the shock front was observed in the simulations, with whistler wave packets 
playing a crucial role. In the first stage of the reformation cycle, a small-amplitude whistler 
perturbation upstream of the ramp is formed This perturbation grows and moves towards 
the ramp. When its amplitude exceeds that of the ramp, this disturbance begins to play the 
role of a new ramp, while the old one moves away downstream. The experimental results 
shown in Figures [30] and [3T] resemble 4 different snapshots for the same shock undergoing 
the reformation. The strongest evidence favouring the shock reformation comes from the 
CIS experiment, which measures the ion composition and full three-dimensional distribu- 
tions for major ions with energies up to 40 keV/e (Rem e et al.|1997| l. The time resolution of 
these measurements is about one spacecraft spin, 4 s. Figure |33] shows 8 snapshots obtained 
at the upstream edge of the shock foot, where the disturbances of the solar wind magnetic 
field are still small. The Figure shows the number of counts vs a function of V x and V y 
in the GSE coordinate system; with the data being integratedin the V z direction. Reflected 
ions are observed for the first time at 07:05:16 (see the maximum of the number of counts 
in the quadrant corresponding to V x < and Vy < in the first snapshot). In the time 
interval from 07:05: 16 to 07:05:44, the position of this maximum in the velocity space does 
not change considerably. In addition, there exists a second population of reflected ions in 
the quadrant corresponding to V x < and V y > 0. From the snapshots it is easily seen 
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Fig. 33 (top) Ion velocity distributions obtained from CIS measurements within the forward part of the foot 
for the Earth's bow shock crossing on 24 January 2001 aboard SCI and (bottom) temporal variations for 
the relative number of counts corresponding to reflected ions. The distributions were calculated in the GSE 
coordinates. In the bottom panel, a blue line corresponds to ions with V y < 0, while a red line shows the data 
for V y > 0. Strong variations of n r , especially for V y < 0, show that the reflection of ions is bursty. The 
relative positions where the measurements were made are indicated by the dots on the magnetic field profiles 
shown as inserts. (Adapted from Lobzin et al. 2007). 



that the numbers of counts corresponding to the reflected ions show approximately periodic 
variations with a very large modulation depth and a period of about 8 s which corresponds 
to half of the proton gyroperiod T^i - To confirm this statement, we performed a summation 
of the number of counts corresponding to these populations, the results are approximately 
proportional to the corresponding number densities, n r . The temporal evolution of these 
number densities normalized with respect to the corresponding maximum values for the 
time interval considered is shown in Figure [33] (bottom). The quasiperiodic variations seem 
to be more pronounced for the first population (blue line), with the minimum-to-maximum 
ratio being as low as ~ 3%. The number of counts for the second population also varies with 
approximately the same period, in phase with that for the first one. It is worth noting that 
the minimum number of counts corresponding to the reflected ions in this region is greater 
than the 'background noise' by a factor of 5, far beyond experimental errors, while for the 
maximum number of counts this factor is as large as 30 if the 'noise' level is estimated in 
the unperturbed solar wind just before the shock crossing. The observed peculiarities of the 
ion dynamics resemble the features found in the numerical simulations of Krasnoselskikh 
|et al.| ( [2Q02| >, where a quasiperiodic front reformation was observed for quasiperpendicular 
shocks with high Mach numbers. In particular, when the leading wave train before the ramp 
attained a large enough amplitude, a new population of reflected ions appeared upstream 
of the precursor. In other words, the reflection of ions is not stationary. It is quasiperiodi- 
cally modified during the reformation process. In this case a spacecraft that moves slowly 
across the shock, will observe the quasiperiodic appearance/disappearance of reflected ions, 
in accordance with experimental results outlined above. 
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7.3 Conclusions 

In this Section we have presented a set of experimental results for a high-Mach-number 
(Mf = 5) quasiperpendicular (#b„ = 81°) bow shock crossing observed by Cluster space- 
craft on 24 January 2001 at 07:0507:09 UT. The structure of this shock gives a clear evi- 
dence of its nonstationary behavior. In particular, the magnetic field profiles measured by 
FGM experiments onboard different spacecraft differ considerably from each other. This 
difference is clearly seen for large-amplitude oscillations, which have relatively short scales 
of about 1-2 s and resemble nonlinear whistler soliton-like structures that is confirmed by 
analysis of their hodograms. WHISPER measurements reveal the presence downshifted os- 
cillations within the electron foreshock, with nonmonotonic variations of their central fre- 
quency, the characteristic time for these variations is comparable with the proton gyroperiod, 
Tgj = 15.5s. From the analysis of data from CIS experiment it follows that the reflection of 
ions from the shock are also highly nonstationary. Moreover, it is shown that the reflection 
is bursty and the characteristic time for this process is also comparable with the ion gyrope- 
riod. From numerous numerical simulations of quasiperpendicular shocks it is well-known 
that for high Mach numbers the shock becomes nonstationary. Moreover, front reformation 
can take place with a characteristic time comparable with the ion gyroperiod. The combi- 
nation of the features outlined above for the bow shock crossing under consideration is the 
first convincing experimental evidence favoring the shock front reformation. 



8 Conclusions 

There exist several models of quasiperpendicular high Mach number shocks. Theoretical 
considerations and computer simulations on todays level are not capable to describe cor- 
rectly all physical process that determine different aspects of shock physics. The only pos- 
sibility to ensure that the theory or modelling correctly capture major physical effects is 
to rely on analysis of experimental data of direct in situ measurements onboard satellites. 
The best adapted for this goal are Cluster satellites since they allow to distinguish spatial 
and temporal variations and during the mission they had different intersatelllite distance that 
allows one to probe the shock on different scales. The difficult task in such investigation 
program consists in formulation of the right questions to be addressed to data and to their 
analysis. Our aim was to determine major physical processes that define characteristics of 
the most important part of the shock front, its ramp and wave activity around it. The prob- 
lems closely related to this major problem are electron heating mechanisms and transition 
of shock behaviour from stationary to nonstationary. We left beyond the scope of our Re- 
view many questions. One can mention ripples, remote sensing of the shock by field aligned 
beam, instabilities behind the shock front. We restricted ourselves by the analysis of scales 
of magnetic and electric fields, the scale of electron heating, determination of the source and 
generation mechanism of precursor whistler wave train and direct observation of the shock 
front reformation by Cluster satellites. Huge collection of data of statistical analysis and of 
the studies of individual events (case studies) leads to the conclusion that the ramp region of 
supercritical quasiperpendicular shock is nothing more than an intrinsic element of nonlin- 
ear dispersive wave structure slightly modified by reflected ions. This interpretation allows 
one to explain in a natural way the whole collection of data that we reported in this Review 
and to understand the transition from stationary to nonstationary shock behaviour when the 
Mach number exceeds nonlinear whistler critical Mach number. The role of anomalous re- 
sistivity is shown to be relatively weak with respect to effects of dispersion. Our study also 
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points out several important opened questions. Presumably the most important is what is 
detailed mechanism of the electron heating and isotropization. Certainly, the evaluation of 
the role of anomalous resistivity can not be considered as the solved issue, the data set used 
is too poor to come to definite conclusions, thus this study still waits new measurements. 
We did not address the problem of particle acceleration and new results presented here will 
certainly have an impact on re-consideration of this important problem. 



A Remark on comparison of computer simulation results with experimental data 

Recently Comi§el et al. ( 201 1 1 made an attempt to perform computer simulations that can properly reproduce 
the realistic physical conditions corresponding to observations. They modelled the shock dynamics using ID 
PIC code with the realistic ion to electron mass ratio under conditions corresponding to shock conditions 
on 24th January 2001 that was observed by Lobzin et al. (2007}. The only difference between the model 
and real plasma parameters is an unrealistic ratio of (u) pe / n ce ). The modelling results clearly showed that 
the shock indeed is nonstationary. However, it was found that there are some important differences between 
the results of the simulations and observations. The major differences can be summarized as follows. The 
electric fields observed in simulations in the close vicinity of the shock front were much higher than the 
electric fields experimentally registered. The energy flux of waves observed in the foot region upstream of 
the shock front was found to be directed toward downstream that clearly indicates that the waves observed 
in simulations are generated by the beam of the reflected ions and not by the ramp region as the dispersive 
mechanism predicts. This gives an indication that the properties of waves observed are much closer to short 
scale lower hybrid or lower hybrid drift waves described in Section|6]and not to those described in Sectons|2] 
and[5] The question arises where does this difference come from. To answer this question one should consider 
some scaling properties of equations describing dynamics of the shock. In order to do that let us re-write our 
equations in dimensionless form making use of natural variables 

v = v/Va, t = tO c i,r = rQ c i/VA, b = B/ Bo, e = qE/ Q c iMj\VA,n = n/no 

The system have several dimensionless parameters that remain and should be taken into account. These are 
rj = (the authors would prefer the letter fi but it is already used for magnetic permeability), \ = ff-, 
and certainly t: ,i • 

To account for the principal difference between the real physical conditions and simulations let us con- 
sider where the parameter \ mav P lav an important role. In dimensionless variables it appears in two Maxwell 
equations 

divli = — — ffij — n e ) 
X 2 

rotB = (riiVi - n e v e ) H 

X 2 8t 

To clarify its role one can consider the properties of linear waves. One can note that the ratio of electric 
to magnetic field is determined by the refractive index of waves. By definition it is 

and it is easy to see that it is exactly this ratio is used in determination of the electric to magnetic field 
ratio. In SI the system in dimensional variables reads 



rot 



' at 



that leads to following estimate for linear waves: 

[t X it] = wit. 

This can be re-written as follows: 

— = - = N. 

E k± ui 
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It can also be expressed in terms of phase velocity 

V ph 

If we take the waves having velocities close to the shock front velocity (approximately standing whistlers 
in a shock front reference frame) the velocity in the plasma reference frame is V up = V sw = MaVa, thus 
the refraction index is 



N 



1 



V 3W MaVa M A n ci M A n ce y/fj 



Be c ( u>„ 2.7 ■ 10 4 

AT cxp = = ( — = s- ~ 230 

P E k± V sw \Q C£ 1.2 ao 2 

where E^^_ is the electric field component perpendicular to the k- vector. 

On the 24th of lanuary 2001 the solar wind velocity was V sw s\y = 440kms . 

7V cxp « 700 

The maximum value of ratio jf- in simulations is 8 thus 

Nsim — 23. 

it is approximately 30 times smaller than in experiment, that means that for the same level of fluctuations of 
the magnetic field the electric field fluctuations are 30 times stronger than in experiment. 

According to our analysis of dimensionless parameters another important difference consists in similar 
overestimate of electric fields due to even small deviations from quasi-neutrallity. One can suppose that 
this can lead to artificial increase of the role of quasi-electrostatic instabilities of short scale lower hybrid 
waves. As a result the dominant waves observed in simulations are similar to those reported in the Section 
"Anomalous resistivity", namely drift lower hybrid type waves. Presumably, the overestimation of the role 
of the electric field and consequently of short scale oscillations and consequently underestimation of the 
role of lower frequency standing precursor whistler waves results in the difference between observations and 
simulation results. If so, the simulated shock is really resistive while the observed one is certainly dispersive. 
To evaluate the influence of this overestimate of the electric field let us evaluate the electric field needed to 
reflect upstream ion flow assuming that for efficient reflection the potential should be of the order of half 
energy of the incident ion. Reflecting potential in nonlinear wave on the scale about 

c R mi V 2 o V 2 

L = 5 = 0.9 ■ 10 6 cm, E ions = — = 0.5 ■ 10 9 ev— = IkeV 

uj pe 2 c z 

This corresponds to the value of the electric field 



0.5 ■ 10 3 V 
2 ■ 0.9 ■ 10 4 n 



30mVm" 



0.7.10^.60V/m.lO^ =i5 iQ . 4 ^ i5nT 
2 • 2 ■ 10 8 

Id fluctuatioi 

ison with the electric field measurements. These effects are illustrated on the Figur j34] where left hand panels 
show electric and magnetic field fluctuations in units similar to those experimentally observed, and right hand 
panels show the data obtained by Cluster satellites for similar parameters (Mach number, angle and P). The 
magnetic field fluctuations in simulations that will be associated with similar electric fields could be 30 times 
smaller. Thus electric fields capable to trap and reflect ions are associated with the magnetic field fluctuations 
that are quite small, namely, less than ~ InT. The ion trapping and reflection can occur in small amplitude 
oscillations in the foot region that can not happen in real shock. Crucial change of ion dynamics certainly 
results in change of the characteristics of the shock front and wave activity around. 

The goal of this remark is not to understate the role and importance of computer simulations for the shock 
studies. We would like to point out that direct comparison of simulation results with the observations needs 
special attention and analysis of the simulation conditions to ensure that the process is properly described. 
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Fig. 34 Comparison of electric and magnetic fields observed on 24th of January 200 1 by Cluster satellites and 
obtained in computer simulation by Comi§el et al. 1 201 1 1. Both panels show the fields in the vicinity of the 
foot/ramp regions. Top panels represent magnetic fields as measured by FGM instrument and STAFF search 
coil magnetometer (right), and obtained as a result of simulations(left). Bottom panels represent measured and 
simulated electric fields. One can see huge difference in amplitude of electric fields between measurements 
and simulations that results from artificial ratio of plasma to gyrofrequency and consequently unrealistic 
refractive index of waves. (Figure is provided by J. Soucek and H. Comisel) 
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B Table of notations used in the article 



Parameter 


Interpretation 


B 


the magnitude of the magnetic field 


Bo 


upstream magnetic field 


B n 


magnetic field component along the normal to the shock 


a 8irnT 
P ~ B i 


the ratio of total particle thermal pressure to the magnetic field pressure 


AB 


change of the magnetic field through inhomogeneous layer 


Mij 


time difference of observation of shock front features such as electric field 


spikes by different satellites i and j 


AE 


an amplitude of the electric field spike feature 




electric field components along corresponding axes 


B spike 


maximum amplitude in electric field spike 
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Parameter 


Interpretation 




electric field as measured by means of probes i and j 


E 


electric field vector 


e ij 


electric potential difference between probes i and j onboard single satellite 


/ = UJ/27T 


wave frequency 


fe = i?ce/27T 


electron gyrofrequency (in Hz) 




growth rate of lower hybrid waves 


k 


wave-vector of a wave 


kij 


an estimate of the k-vector component from electric field probes i and j 




measurements onboard one single satellite 


fe,. 


parallel to the magnetic field component of the wave- vector 


Igr 


characteristic gradient scale inside the inhomogeneous layer 


L d 


dissipative scale 




dispersive scale 




ion, electron inertial length 


LBr 


thickness of the ramp region of the shock as seen in magnetic field mea- 




surements 




characteristic scale of the electrostatic potential variation in the shock front 


Lan 


characteristic scale of energy exchange due to anomalous resistivity 


Lf 


characteristic spatial size of the magnetic foot 


L r 


width of the magnetic foot 




ion, electron mass 


M A = Vup/V A 


Alfvenic Mach number, the ratio of the normal component of the upstream 




flow velocity to Alfven velocity 


Mms = Vup/Vms 


magnetosonic Mach number 


M w = Vw,max /Va 


nonlinear critical whistler Mach number 


n E 


shock front normal determined from timing of electric field spikes mea- 




sured onboard four satellites 


n B 


shock front normal vector from magnetic field measurements 


n 


plasma density 


n 


shock front normal vector 


n r 


number density of reflected ions 




convective ion gyroradius 




electron Larmor radius 


Sy 


Poynting flux along the magnetic field 


T 


total plasma temperature 


T 


ion, electron temperature 


T e \\ 


parallel to magnetic field electron temperature 


T el _ 


electron temperature perpendicular to the magnetic field 




ion gyroperiod 




the relative velocity of electrons carrying current 


V A 


Alfven velocity 




the normal component of the upstream velocity to the shock surface in its 




rest frame 


V Ms 


velocity of the magnetosonic wave propagating in the same direction as 




the shock to the background magnetic field 




highest possible velocity of nonlinear whistler wave that can stay in the 




upstream flow 


Vsw 


solar wind velocity 


Vss 


relative shock spacecraft velocity 


Vph 


phase velocity of ion sound wave 


v ph 


phase velocity of wave 


Vis 


ion sound velocity 


Vx,y,z 


velocity components along corresponding axes 


w k 


electric field energy density 


X 


wavelength of precursor whistler wave 


Ad = e k B T e /ne 2 


Debye radius 




effective collision frequency due to wave particle interaction 


6bu 


the angle between the magnetic field and shock front normal 


OkB 


angle between the magnetic field vector and the wave vector 
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Parameter 


Interpretation 


^Ih ~ V ^ci^ce 

Q ci , e = eB/mi ;e 

OJp 


lower hybrid frequency 
ion plasma frequency 
ion, electron gyrofrequency 
electron plasma frequency 
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